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ABSTRACT 

Laser  produced  plasma  jets  interacting  with  a  background  plasma  have 
been  used  to  study  magnetic  field  generation  in  shock  waves.  Shock  heating 
provides  axial  electron  temperature  and  density  gradients  which  are 
perpendicular  to  the  radial  temperature  and  density  gradients.  Electron 
heat  transport  and  ion  diffusion,  in  the  radial  direction,  occur  at  different 
rates.  This  combination  of  non-parallel  temperature  and  density  gradients 
generates  a  magnetic  field  in  the  azimuthal  direction.  Simulations 
corroborate  the  experimental  observation  that  magnetic  fields  are  generated 
when  a  supersonic  plasma  jet  interacts  with  a  background  plasma.  Magnetic 
flux  generated  by  this  mechanism  requires  no  initial  field,  which  is  in 
contrast  to  the  dynamo  mechanism  which  requires  an  initial  seed  field. 
Specific  applications  analyzed  in  this  thesis  include  interplanetary  shocks 
and  nuclear  EMP  effects  in  the  MHD  domain.  It  may  be  assumed  that 
shock  generated  magnetic  fields  are  also  of  importance  under  astrophysical 
conditions. 
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I.  INTRODUCTION 

Since  the  early  1960’s,  extensive  research  has  been  done  on  laser  pro¬ 
duced  plasmas.  The  motivation  for  this  stemmed  from  the  fact  that  laser 
produced  plasmas  provide  an  efficient  and  cost  effective  way  of  studying  col- 
lisional  plasma  phenomena.  Fusion  research  has  drawn  greatly  from  the 
knowledge  obtained  from  laser  plasma  research.  Inertial  and  magnetic 
confinement  fusion  techniques  are  prime  examples  of  this.  Particularly  in 
the  field  of  inertial  confinement  fusion,  a  thorough  understanding  of  how 
intense  coherent  laser  light  interacts  with  the  surface  of  a  deuterium  pellet 
is  critical  in  achieving  the  enormous  temperatures  and  densities  required  to 
obtain  an  appreciable  fusion  yield.  Laser  plasma  surface  interactions  have 
been  extensively  studied  at  the  Naval  Postgraduate  School. 

High  density  plasma  jets  are  produced  when  the  incident  laser 
radiation  strikes  an  aluminum  target.  When  the  laser  light  initially  hits  the 
target,  electrons  within  the  aluminum  absorb  the  incident  radiation.  Since 
the  collision  frequency  of  electrons  in  the  target  is  on  the  order  of  1013  sec-1, 
within  the  laser  pulse  length  of  22  nanoseconds,  the  electrons  have  more 
than  ample  time  to  distribute  their  energy  to  other  electrons  as  well  as 
lattice  phonons,  hence  allowing  the  optical  lase?  energy  to  be  converted  to 
thermal  energy.  The  highly  ionized  target  material  forms  a  dense  plasma 
jet  which  expands  in  the  direction  that  is  normal  to  the  target  surface. 
During  heating  and  ionization  of  target  material,  radial  temperature 
gradients  associated  with  the  laser  radiation  and  axial  density  gradients, 
which  are  associated  with  the  expanding  dense  target  material,  generate 
magnetic  fields  which  have  the  structure  shown  in  Figure  1.  If  the  plasma 
jet  is  viewed  head  on,  the  laser  produced  magnetic  field  direction  is 
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clockwise.  These  fields  decay  relatively  slowly  when  compared  to  the  time 
scale  of  jet  expansion  into  the  background  plasma.  Jet  expansion  occurs  in 
<  1  |isec  over  a  distance  of  3  cm.  During  this  expansion  the  magnetic  fields 
are  "frozen  in"  the  expanding  plasma  jet.  The  magnetic  fields  have 
azimuthal  symmetry  with  respect  to  the  axis  of  the  jet.  In  this  geometry 
the  magnetic  field  lines  closely  represent  concentric  rings  around  the  axis  of 
the  expanding  plasma  jet. 


Figure  1:  Geometry  of  Laser  Produced  Plasma  Jet  and  Associated  Mag¬ 
netic  Field. 

The  fields  and  the  motion  of  ions,  in  direction  normal  to  the  target 
surface,  tend  to  prevent  expansion  in  the  radial  direction;  hence  the  ionized 
target  material  closely  represents  a  plasma  jet  expanding  outward  at 
approximately  106  m/sec.  Comparison  of  this  velocity  with  that  of  the  ion 
acoustic  velocity  of  the  background  plasma  immediately  indicates  that  the 
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jet  is  supersonic.  Calculations  will  be  shown  in  this  thesis. 

Once  the  laser  pulse  has  ceased,  additional  magnetic  fields  are 
generated  by  nonparallel  temperature  and  density  gradients  created  by  a 
shock  transition  at  the  leading  edge  of  the  jet.  If  the  background  plasma 
density  is  sufficiently  high,  a  shock  front  is  created  at  the  leading  edge  of 
the  plasma  jet  which  causes  a  discontinuity  in  the  electron  and  ion 
temperatures  in  accordance  with  classical  shock  theory.  Moreover,  since  the 
jet  is  highly  supersonic,  temperature  is  increased  by  an  order  of  magnitude 
higher  than  density.  This  is  because  at  higher  densities  the  plasma  (as  with 
any  gas)  tends  to  become  more  difficult  to  compress  (see  Figure  2).  These 
shock-generated  magnetic  fields  have  the  geometry  as  shown  in  Figure  3. 
Notice  that  the  field  direction  is  counterclockwise  when  viewed  head  on. 


Mach  Number 

Figure  2:  Temperature  and  Pressure  Characteristics  Versus  Mach 
Number  assuming  no  ionization  losses. 
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Figure  3:  Geometry  of  Shock  Produced  Magnetic  Fields. 

Once  the  plasma  has  been  shock  heated,  the  hot  magnetized  electrons 
(due  to  their  small  Larmor  radius)  rapidly  disperse  their  thermal  energy 
through  electron-electron  collisions.  The  ions  have  a  relatively  large  Larmor 
radius  and  mean  free  path  ,  rL(f» 8^,  8^  is  the  characteristic  radial 

size  of  the  jet).  The  hot  ions,  therefore,  expand  out  radially  at  their 
respective  thermal  (sound)  velocity.  Diffusion  of  electron  thermal  energy  via 
electron-electron  interactions  (while  in  their  gyromagnetic  orbits)  occurs 
faster  than  the  radial  expansion  of  the  ions.  Therefore,  at  later  times  the 
radial  electron  temperature  gradient  has  a  characteristic  width  (half  width 
at  half  maximum)  which  is  greater  than  the  radial  electron  density  gradient. 
The  electron  density  gradients  both  in  the  axial  and  radial  directions  are 
inherently  dictated  by  the  ion  density.  Hence  the  radial  electron  density 
gradient  can  not  relax  as  fast  as  the  radial  electron  temperature  gradient. 
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That  is,  the  radial  electron  density  distribution  in  the  laser  plasma  does  not 
change  as  rapidly  as  the  electron  temperature.  Electrons  colliding  with 
other  electrons  while  undergoing  gyromagnetic  motion  in  the  laser 
generated  magnetic  field,  can  transmit  energy  much  more  quickly  than  ions 
can  provide  particle  transport.  This  is  mentioned  here  as  point  of  fact, 
calculations  will  however  be  shown  later  in  this  thesis. 

The  study  of  these  dense  plasma  jets  and  their  surrounding  magnetic 
fields  is  important  not  only  in  perfecting  inertial  confinement  fusion,  but 
also  in  understanding  the  magnetic  field  structure  of  plasma  jets  which  are 
directly  related  to  phenomena  such  as:  electromagnetic  pulse  effects 
associated  with  nuclear  detonations,  magnetic  fields  associated  with  the 
earth’s  bow  shock,  matter  jets  associated  with  some  galaxies  and  quasars, 
and  other  occasions  where  plasma  jets  exist  in  astrophysical  phenomena. 
Chapter  VI  of  this  thesis  will  provide  some  additional  study  in  plasma  jet 
related  phenomena. 

A.  THESIS  STATEMENT 

Shock  heating  at  the  front  of  a  plasma  jet  creates  temperature  and 
density  gradients  which  generate  unique  magnetic  fields.  They  are  unique  in 
that  their  direction  is  reversed  when  compared  to  the  self-generated 
magnetic  fields  which  are  associated  with  laser-surface  interactions. 
Numerical  simulations  have  been  used  to  show  that  shock  heating, 
combined  with  fast  thermal  transport  in  the  radial  direction,  generates  the 
experimentally  observed  shock  magnetic  fields.  Analytical  work  with 
applications  to  interplanetary  shocks  and  nuclear  EMP  in  the  MHD  domain 
is  also  presented. 
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II.  THEORY 


A.  SELF-GENERATED  MAGNETIC  FIELDS 

The  magnetic  fields  previously  discussed  are  inherently  self-generated. 
That  is,  they  require  no  initial  seed  magnetic  field.  For  example,  the 
terrestrial  magnetic  field  exists  today  due  to  a  process  called  "the  dynamo 
mechanism".  This  mechanism  requires  an  initial  seed  magnetic  field. 
Nonuniform  rotations  within  the  earth’s  molten  core  (driven  by  irregular 
convection  of  conducting  fluid)  greatly  slow  the  decay  of  the  initially  present 
magnetic  fields,  while  continuing  to  generate  additional  magnetic  field  flux 
by  ifxljt  induced  currents.  The  theory  of  magnetic  field  generation,  which 
will  be  described  here,  is  strikingly  different  in  this  respect.  The  important 
point  about  the  theory  presented  here,  is  that  a  source  term  exists  in  the 
magnetic  field  generation  equation  which  allows  spontaneous  growth  of 
magnetic  flux  in  a  plasma. 

In  order  to  generate  a  magnetic  field  it  is  immediately  apparent  that 
there  must  exist  an  electric  field  which  satisfies  Maxwell’s  equation: 

|=-Vx^.  (1) 

at 

Generalized  Ohm’s  law  provides  an  equation  of  motion  for  the  electrons 
which  respond  to  the  influence  of  an  electric  field.  Generalized  Ohm’s  law  is 
given  by: 

-  *■ 

pe  ^-  +  (l?e-V)l?e  =  -ene(£  +  D>ex#)-VPe  +  -^J\  (2) 

at  a 

The  subscript  e  denotes  electron  parameters.  pe  indicates  mass  density,  j 
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indicates  current  density,  a  indicates  electrical  conductivity,  and  if  indicates 
fluid  velocity.  The  electron  plasma  frequency: 


e>p  = 


nee‘ 


torn 


v* 


or  vp  =8.99  We*  where  vp  is  given  in  sec-1  and  ne  is  in  m-3,  is  about 
l.lxlO14  sec-1.  Therefore,  it  is  justifiable  to  neglect  the  inertial  and  convective 
terms  on  the  left  hand  side  of  equation  (2).  Thus  equation  (2)  becomes, 


g  =  lj-0fe  x i?) - — VPe  . 


n.e 


(3) 


Taking  the  curl  of  equation  (3)  and  using  Maxwell’s  equation  given  in 
equation  (1)  gives, 


Mz-IvxJ+Vx^x^  +  Vx  TTTVP«  • 


nee 


(4) 


Now  by  using, 

H0/  =  V  x  £?  and  V  •  2?  =  0 

the  first  term  on  the  right  hand  side  of  equation  (4)  becomes, 

_  i-  V  x/=  —  V2#  . 
o  \loO 

The  equation  of  state,  assuming  an  ideal  gas,  Pe  =  nekTe  gives, 


(5) 


(6) 


V  x  — VPe  =  V  x  -^-V(nekTe)  =  V  x 


nee 


nte 


—  VTe  +  —  Vne 
e  n.e 


(7) 


Since  V  x  VTe  =  0,  equation  (7)  may  be  rewritten  as, 


V  x  — VPe  =  V  x 


nee 


kTe 


Vn. 


nee 


(8) 


Now  using  the  vector  identity, 

V  x  (<j^)  =  $V  x  X  +  V<|>  x  X, 
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with  the  scalar  function  <)>  =  kTe!nee,  the  vector  X  =  Vne ,  and  V  x  Vne  =  0,  the 
right  hand  side  of  equation  (8)  becomes, 

=  *-VTe  x  -1-V/v  (9) 

e  ne 

Finally  equation  (4)  becomes, 

-^  =  —  V2#  +  Vx(U»c  x&)+-VTe  x  —  Vne.  (10) 

dt  jx0o  e  ne 

The  first  term  on  the  right  side  of  equation  (10)  represents  diffusion  of  the 
magnetic  field  through  the  plasma.  The  rate  of  field  diffusion  or  field  decay 
in  the  plasma  will  be  determined  by  the  electrical  conductivity,  or  Spitzer 
resistivity,  of  the  plasma.  Based  on  dimensional  analysis  it  is  obvious  that 
the  characteristic  time  required  for  magnetic  field  decay  is,  xm  =  |i0a/2  (l 
represents  the  characteristic  size  of  the  magnetic  field  region).  If  Ohmic 
heating  losses  are  negligible,  then  the  current  (i.e.  the  magnetic  field)  will 
not  decay,  and  the  diffusion  term  can  be  neglected.  The  second  term 
represents  the  convective  element  of  changes  of  S  within  the  plasma. 
Equation  (10)  without  the  third  term  is  the  called  "the  dynamo  equation" 
which  when  coupled  with  the  equations  which  describe  the  convective 
behavior  (if)  of  fluids  within  a  rotating  body,  can  be  solved  (often 
numerically)  to  provide  an  idea  of  how  dipoles  and  higher  multipoles  are 
generated  within  the  earth  and  other  celestial  bodies.  The  ratio  of  the 
convection  term  to  the  diffusion  term  is  often  called  the  magnetic  Reynolds 
number,  given  by, 

R  =  \ioOvl  =  -j- , 

where  l  is  the  characteristic  size  of  the  field  or  convective  cell.  For 
experiments  analyzed  in  this  thesis,  the  magnetic  Reynolds  number  is  large, 
implying  that  the  field  is  "frozen"  in  the  plasma. 


V  x 


kT ' 


en. 


Vn, 
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The  third  term  represents  a  source  term  which  allows  for  spontaneous 
magnetic  field  generation  when  non-parallel  temperature  and  density 
gradients  are  present.  It  is  this  term  which  allows  for  the  growth  of  a 
magnetic  field  in  the  absence  of  an  initial  seed  field.  The  initial  laser 
induced  self-generated  magnetic  fields  are  produced  in  this  way.  Laser 
irradiation  of  the  target  surface  rapidly  raises  the  material  temperature  at 
the  surface.  A  radial  temperature  gradient  is  established  due  to  the  energy 
density  cross  section  of  the  laser  beam.  The  expanding,  ionized  target 
material  forms  a  dominant  axial  density  gradient,  which  propagates  at  the 
jet  velocity  away  from  the  target  surface.  The  dominant  radial  temperature 
and  axial  density  gradients  generate  the  initially  observed  magnetic  fields 
as  shown  in  Figure  1.  Subsequent  expansion  of  the  jet  at  substantial 
background  pressures  allows  shock  formation  to  occur,  which  itself  is  a 
mechanism  which  generates  temperature  and  density  gradients,  and  allows 
for  magnetic  field  production.  Shock  formation  provides  heating  and 
compression  behind  the  shock  front. 


B.  SHOCK-GENERATED  MAGNETIC  FIELDS 

In  the  case  of  shock  generated  fields,  differences  in  the  radial 
temperature  and  density  gradients  exist  because  the  radial  electron  density 
gradients  are  maintained  by  the  relatively  slower  moving  ions.  This  must  be 
true  in  order  to  maintain  quasineutrality  of  the  plasma.  Electrons  maintain 
a  force  balance  between  the  attractive  coulomb  interaction  with  the  ions  and 
the  pressure  gradient  force.  Equilibrium  between  these  two  opponents 
ensures  that  quasineutrality  on  scales  larger  than  a  Debye  length  is  obeyed. 
The  Debye  length  is  given  by, 


%kTe 


nee‘ 


Vt 
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or, 


XD  =  740 


Vi 

J 


where  Te  is  in  eV,  ne  is  in  cm-3,  and  XD  is  in  cm.  For  the  plasmas  considered 
here  ne  =  1014  cm-3,  kTe  =  42  eV,  and  hence  Xd  ~  4.8X10-4  cm.  This  is  much 
smaller  than  the  dimensions  of  the  plasma  jet.  Therefore,  charge  neutrality 
on  scales  of  the  jet  size  are  maintained. 

While  the  electron  density  gradients  are  largely  dictated  by  the  slower 
ions,  electron  temperature  gradients  are  not  only  dictated  by  electron 
diffusion  but  also  by  the  heat  conductivity  of  the  magnetized  electrons. 
Conduction  of  heat  by  the  electrons  in  and  around  the  plasma  jet  occurs 
relatively  quickly  compared  to  the  ion  motion. 

Two  distinct  cases  exist  when  studying  plasma  jets.  First,  the  plasma 
jet  may  simply  propagate  into  a  vacuum.  In  this  case,  free  expansion  occurs 
and  from  the  onset  of  plasma  production  the  relevant  gradients 
monotonically  decay,  hence  field  generation  due  to  non-parallel  gradients  is 
limited  and  quickly  becomes  negligible.  Second,  if  the  plasma  jet  expands 
into  a  relatively  dense  background  plasma,  interaction  with  the  background 
plasma  produces  higher  laser  produced  and  shock  generated  magnetic  field 
strengths  at  the  front  of  the  plasma  jet.  Shock  generated  magnetic  fields  are 
identified  experimentally  by  reversal  of  the  azimuthal  field  component  at  the 
jet  front  (see  Figure  4). 

Larger  magnetic  fields  might  be  expected  since  an  ionized  background 
provides  a  highly  conductive  medium  so  that  larger  currents  are  produced 
by  the  existing  gradients.  Shock  generated  fields  occur  when  the  jet 
interacts  with  the  background.  As  the  plasma  jet  propagates  into  the 
background  plasma  at  a  velocity  greater  than  the  local  ion  acoustic  velocity, 
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Figure  4:  Magnetic  Field  Strength  Versus  Background  Plasma  Pres¬ 
sure.  Maximum  azimuthal  magnetic  field  at  z=0.4  cm,  r=0.3  cm,  0=0°, 
as  a  function  of  N  2  background  pressure  for  an  incident  laser  power  of 
300  MW.  [Ref.  1] 

a  supersonic  shock  front  is  created  at  the  leading  edge  of  the  jet.  The  shock 
front  is  a  thin  (thickness  =  3  mm),  nearly  discontinuous,  transition  in 
plasma  properties. 
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C.  SHOCK  THEORY 


The  ion  acoustic  speed  in  a  plasma  is  easily  derived  by  linearizing  the 
ion  equation  of  motion  and  continuity  equation.  These  are  given  by, 


f+w 


V)u> 


=  -VP 


(11) 


and 

It  +  v  •  (pif)  =  o ,  (12) 

ot 

respectively.  In  order  to  linearize  these  equations,  consider  a  uniform  fluid 
with  density  p0,  pressure  P0,  and  velocity  if0  with  small  perturbations  pi,  Pi, 
and  o\,  so  that  p0  +  Pi>  P0  +  Pi,  and  v*0  +  v*i,  describe  the  full  dynamic 
behavior.  Using  p  =  p0+p!,  P  =  P0  +  P1(  and  if  =  \f0  +  in  the  above 
equations  and  taking  Vp0  =  VP0  =  V  •  tf0  =  0  gives, 


(Po  +  Pi)“j^  +  <Po  +  PiXt^i  *  V)i? i  =  -VP  i 


and 


dpi 

dt 


+  0*i  •  Vpx  +  (po  +  Pi)V  •  if  1  =  0. 


Neglecting  second  order  approximations,  i.e.  terms  such  as  piO*x,  provides 
the  linearized,  first  order  approximations  to  the  dynamical  equations: 


ao% 


ap 


ap 

V  J 


Vpx  and 


dpi 

dt 


+  p0V-  0*i  =  0. 


Taking  the  divergence  of  the  first  equation  and  subtracting  it  from  the  time 
derivative  of  the  second  equation  gives, 

d2Pi 


ap 


ap 


V^pi. 


(13) 


Equation  (11)  clearly  represents  the  wave  equation  for  px,  hence  the  ion 
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acoustic  sound  speed  is  given  by, 

5P 

jpJ; 

The  subscript  s  indicates  the  partial  derivative  is  taken  at  constant  entropy, 
and  hence  represents  the  adiabatic  ion  sound  speed  for  small  disturbances. 
The  important  concept  to  take  from  this  presentation  is  that  if  the 
disturbances  are  no  longer  small,  then  second  order  perturbations  can  no 
longer  be  neglected.  Linear  analysis  (as  given  above)  is  no  longer  valid, 
because  the  nonlinear  terms  such  as  pxV  •  V*i  can  not  be  neglected.  Hence, 
shock  wave  formation  is  inherently  a  nonlinear  process. 

Without  exactly  solving  the  nonlinear  fluid  equations,  important 
characteristics  concerning  shock  fronts  can  still  be  derived.  Shock  fronts  in 
fluids  constitute  near  discontinuities  in  T,  p,  and  P.  When  viewed  from  the 
shock  front  frame  of  reference  the  conditions  which  relate  parameters  such 
as  p,  P,T  across  the  shock  front  can  be  more  easily  understood.  In  this 
frame,  the  shock  front  is  stationary  and  the  unshocked  fluid  (denoted  by 
region  1)  is  moving  at  velocity  v  i  while  the  shocked  fluid  (denoted  by  region 
2)  is  moving  at  velocity  v2-  Therefore,  unshocked  and  shocked  quantities  will 
be  denoted  with  subscripts  of  1  and  2,  respectively.  Figure  5  shows  this 
arrangement. 

Three  quantities  must  be  conserved  when  passing  through  the  shock 
front:  (1)  mass  flux,  (2)  momentum  flux,  and  (3)  energy  flux.  This  assumes 
that  the  shock  triggers  no  local  energy  release  (such  as  by  ionization).  As  it 
turns  out,  this  assumption  is  not  valid  for  the  strong  shocks  considered  in 
this  thesis.  The  presence  of  ionizations  will  be  discussed  in  subsequent 
sections  in  order  to  explain  discrepancies  between  observed  and  simulated 
magnetic  field  magnitudes.  First,  mass  conservation  across  the  shock  front 
follows  directly  from  the  continuity  equation.  If  the  continuity  equation  is 
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T  /  p  /  V  ,  c 
2  2  2  2 

< - 


Shock  Front 

Figure  5:  Shock  Front  Reference  Frame. 

integrated  over  the  shock  front  then, 

u  lP  1  =  u  2P2  =  J  »  (14) 

where  vt  p£  represent  the  velocity  and  density  in  the  region  i ,  respectively.  J 
is  the  mass  flow  per  unit  area  per  unit  time  across  the  shock  front.  This 
equation  states  simply  that  mass  can  neither  be  created  nor  destroyed  at 
the  shock  front.  Second,  if  the  continuity  equation  and  equation  of  motion 
are  combined  and  integrated  over  the  shock  front,  then  the  momentum  flux 
density  for  the  flow  normal  the  shock  is  conserved,  that  is: 

Pi  + Pivf  =P2  + P2V2  >  (15) 

where  the  term  pu2  is  the  momentum  flux  density  due  to  a  fluid  element, 
and  the  pressure  term  represents  the  momentum  flux  associated  with  ion 
and  electron  thermal  motion.  Third,  the  energy  flux  into  the  shock  must  be 
the  same  as  the  energy  flux  flowing  out  of  the  shock.  Mathematically,  this 
is  equivalent  to  Bernoulli’s  principle  which  states  that  for  compressible 
fluids: 


T  ,  p  ,  V  ,  c 

1111 

< - 
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(16) 


1  «  Pi  U  i  1  o  P2 

2  Pi  Pi  2  p2 


U  2 

P2  ’ 


where  u*  represents  the  specific  internal  energy  of  the  fluid  in  region  i 
(t  =1,2).  Only  those  fluid  elements  with  components  of  velocity  normal  to  the 
shock  front  undergo  a  discontinuous  jump  in  density,  temperature,  and 
pressure.  It  is  assumed  at  this  stage  that  any  magnetic  fields  which 
permeate  the  plasma  jet  do  not  affect  the  velocity  components  which  are 
tangential  to  the  jet  or  shock  front.  Hence,  the  tangential  velocity 
components  remain  continuous  across  the  shock  front.  By  introducing  the 
enthalpy  H  =  U  +  PV,  it  is  possible  to  solve  for  an  expression  which  only 
depends  on  the  initial  and  final  thermodynamic  state  of  the  plasma  and  is 
independent  of  the  jet  velocity,  v  i.  Using  equation  (16)  and  h  =  u  +P,  where 

h  =  y  gives, 


Next  using  mass,  momentum,  and  energy  conservation  to  find  an  expression 
independent  of  v  x  and  v  2  gives, 


hi 

Pi 


1 

2 


Pi  +  P2 
PlP2 


(P2~P!)  =  0. 


Or  using  the  explicit  quantities  gives, 

H1-H2=  |(P1  -  P2XVx  +  V2).  (17) 

Equation  (17)  defines  the  shock  adiabat,  or  the  so  called  "Hugoniot".  Notice 
that  this  expression  is  independent  of  all  dynamic  variables,  and  hence 
defines  a  static  relationship  between  the  shocked  and  unshocked  properties 
of  the  plasma.  From  examination  of  equation  (17),  it  is  seen  that  since 
P2  >P  1,  the  shocked  plasma  enthalpy  must  be  greater  than  the  unshocked 
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enthalpy.  The  shock  front  of  a  plasma  jet  with  cylindrical  symmetry  acts  to 
create  temperature,  density,  and  pressure  gradients  in  the  •£  and  -t 
directions.  If  Mx  :>  2.4,  then  the  shock  front  raises  temperature  much  more 
efficiently  than  density,  as  seen  in  Figure  2.  Therefore,  the  generated 
temperature  gradient  will  become  larger  than  the  density  gradient  in  the  £ 
direction. 

Equations  (14),  (15),  (16)  and  the  perfect  gas  law, 

h  =c„T  =  — —  where  y=  — , 
p  Y- 1  P  c„ 


can  now  be  used  to  derive  all  the  classical  relationships  between  the 
upstream  and  downstream  properties  of  the  plasma.  For  example: 


p2  =  (Y+lWi 
Pi  (y  -  1)M  x  +  2  ’ 

P2  1 +  yMi 
P  i  l+Y^f’ 


~  =  [2  -Ky-lW?] 
1 1 


2yAf  i  -  (y-  1) 
(y  +  1)2M  i  ’ 


(18) 

(19) 

(20) 


where  Af1  =  u1/csl  and  M2  =  v2/cS2  are  the  upstream  and  downstream  Mach 
numbers  respectively.  Equations  (18)  and  (20)  are  plotted  in  Figure  2  for 
Y  =  1.66  which  corresponds  to  a  monatomic  gas. 


Experimentally,  the  background  nitrogen  gas  becomes  ionized  by  the 
radiation  emitted  from  the  hot  laser  produced  plasma  (=  100  eV).  The 
electron  temperature  of  the  ionized  background  is  assumed  to  be  on  the 
order  of  1.0  eV.  Ion  acoustic  velocity  is  calculated  by: 

IV4 


C,  = 


yZkTe 


mi 


=  9.79x10s  eV 


(21) 


where  c8  is  the  ion  acoustic  velocity  in  cm/sec,  p  is  the  ratio  of  ion  to  proton 
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mass,  Z  is  the  ionization  state  of  the  gas,  and  Te  is  the  electron 
temperature  in  eV.  Bird  [Ref.  2]  experimentally  measured  the  jet  velocity  to 
be  approximately  2.8xl06  cm/sec,  with  experimental  error  of  about  ±10%.  The 
plasma  jet  Mach  number  is  a  function  of  the  nitrogen  background 
temperature.  Using  y=1.66,  Z- 1,  and  |a.=28  gives: 


Mx 


11.75  eV* 

VTI 


where  Tx  is  the  background  temperature  in  eV.  From  equation  (20)  it  is  also 

rp 

evident  that  the  Mach  number  is  a  function  of  — ^ -.  Hence,  T2  is  a  function 

T  i 

of  Tx.  This  functional  dependence  is  given  by: 

T2  =  -1.35xl0~3  eV1  T?  +  0M6TX  +  42.75  eV  . 


A  background  temperature  of  1.0  eV  then  implies  that  the  shock  heated 
temperature  is  43.6  eV,  and  the  Mach  number  (Mx)  is  11.75.  This  is  a  fairly 
strong  shock,  which  implies  that  further  excitations  and  ionizations  will 
occur.  N2  and  N2  have  dissociation  potentials  of  9.758  eV  and  8.72  eV, 
respectively.  Diatomic  nitrogen’s  first  ionization  potential  is  15.58  eV. 
Monatomic  nitrogen  has  ionization  potentials  of  14.534  eV  and  29.601  eV 
[Ref.3].  Likewise,  aluminum  has  ionization  potentials  of  5.986  eV,  18.826 
eV.  and  28.448  eV.  All  these  atomic  processes  are  possible  mechanisms  of 
energy  loss  from  the  shock  heated  plasma.  Energy  used  for  excitations  and 
ionizations  will  not  contribute  to  raising  electron  temperatures  to  the 
theoretical  level.  Reduced  thermal  gradients  imply  smaller  magnetic  fields. 
The  initial  analysis  completed  here  will  neglect  these  losses  and  hence  will 
yield  magnetic  field  magnitudes  which  are  too  high.  Assuming  that  shock 
heating  of  the  electrons  is  on  the  order  of  10%  efficient,  magnetic  field 
magnitudes  do  fall  in  line  with  experiment. 
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The  measured  jet  velocity  given  above  was  obtained  at  a  background 
pressure  of  700  mtorr,  which  was  sufficiently  high  to  cause  shock  formation. 
For  pressures  below  about  200  mtorr,  shock  formation  was  not  observed.  At 
lower  background  pressures  the  plasma  jet  expanded  with  greater  velocity 
(up  to  =  10  times  greater);  however,  dynamic  gradients  rapidly  fall  off  due  to 
free  expansion  so  that  significant  self-generated  magnetic  fields  are  minimal 
and  shock  heating  is  nonexistent. 


D.  ELECTRON  HEAT  TRANSPORT 


Cross  field  electron  heat  conduction  is  a  much  more  complicated 
process  which  is  still  being  studied  today  in  connection  with  fusion 
confinement  applications.  Electron  heat  conduction  across  a  magnetic  field  is 
a  topic  of  continuing  research,  and  one  which  is  not  very  well  understood. 
The  third  moment  of  the  collisional  Boltzmann  equation  describes  the 
behavior  of  heat  flow  in  a  plasma  with  a  3-dimensional  distribution  function, 
f  ),  given  by: 


K 

dt 


v>-vf  +  —  ■  y~  = 

m  du* 


where  P  represents  body  forces  such  as  gravity  or  the  Lorentz  force. 
Multiplying  the  above  equation  by  Y-mtf  and  integrating  over  velocity  space 
yields  the  third  moment  of  the  Boltzmann  equation.  Assuming  the  body 
forces  are  caused  by  a  Lorentz  force  then  the  above  equation  becomes 
[Ref.5]: 


\m  J  +\ m  !  v*(lf  ‘  vf  W  "  f  J  =  -y  j  v2 


dv> 


K 

dt 


dtf. 


Performing  the  integration  yields: 


d 

•  ^ 
nmv2  .  3 nkT 

_  .  n* .  v  •  ft _ sL 

•  * 

nmv 2 

dt 

2  +  2 
k  . 

—  nr  u  t  v  ±i  —  *^7" 
at 

2 

»  M 
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where  it  is  energy  flux  and  is  given  by, 

it  =  -j-fim  J  if  (if  •  i?)f(?,i?,t)di f, 

JU 

and  n  is  the  density  given  by, 

n  =  J  fif(~r*,T?,t)  <My 

where  A  is  the  average  particle  density  given  by  The  heat  flux  depends 

on  the  3-dimensional  distribution  function  of  the  electrons,  which  itself  is 
not  exactly  known  and  must  also  be  solved.  Kinetic  theory  in  this  case  does 
not  provide  a  very  eloquent  solution  to  the  electron  heat  transport  problem. 
The  problem,  however,  becomes  more  tractable  if  the  specific  nature  of  the 
electron  and  ion  motion  is  analyzed.  The  electrons  will  be  examined  in  this 
section.  Ions  will  be  examined  in  the  next  section.  For  this  analysis  it  is 
assumed  that  experimentally  measured  density  values  are  accurate  and 
that  the  plasma  is  100%  ionized.  Therefore,  no  neutrals  are  present.  It  is 
also  assumed  that  ion  density  is  equal  to  electron  density  and  electron 
temperature  is  larger  or  equal  to  ion  temperature  in  any  given  region. 
Although  the  second  assumption  may  not  be  strictly  valid,  deviations  within 
one  order  of  magnitude  do  not  change  the  physics  of  the  mechanism. 
Important  electron  parameters  such  as  mean  free  path  (>.e),  Larmor  radius 
(rLje),  and  the  ratio  of  gyrofrequency  to  collision  frequency  (cocTle)  are 
calculated  as  a  function  of  range  from  the  plasma  jet  axis.  All  variables  are 
in  cgs  units.  Temperature  is  in  eV.  The  electron  mean  free  path  (Xe )  is  given 
by  [Ref.4]: 

K  =  3.4xl013-^^  cm  .  (22) 

n  lnA 

In  A  is  called  the  Coulomb  logarithm  and  is  included  to  account  for  the 
cumulative  effect  of  small  angle  collisions,  since  these  collisions  occur  much 
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more  frequently  than  large  angle  collisions.  A  is  defined  as  the  average 
ratio  of  Debye  length,  \D ,  to  the  impact  parameter  radius,  rQ : 

X 

A  =  — —  =  12: in  Xn. 

r0 


In  laboratory  plasmas,  for  example,  if  n  ~  1011  cm-3,  and  kTe  =2  eV,  then 
In  A  =  10.2.  The  logarithm  of  A  is  relatively  insensitive  to  changes  in  plasma 
parameters.  Assuming  a  gaussian  temperature  and  density  profile  and 
In  A  =  10  gives: 


K 


1.07 


((Tmax-T6g)e"r2/8^ 

max  ^bg  )® 


+  n»>2 

+  ^Ig) 


cm  . 


Tm ax  and  Tbg  are  the  maximum  electron  temperature  and  background 
electron  temperature  in  eV,  respectively.  nmax  and  nbg  are  the  maximum 
electron  density  and  background  electron  density  in  relative  units, 
respectively.  is  the  radial  electron  density  width  and  6^  is  the  radial 
electron  temperature  width.  Now  letting  71max=43.597  eV,  Tbg= 1.0  eV, 
nmav=14.9  r.u.,  5r2r=0.3  cm2,  and  8%  =0.3  cm2  provides  the  curve  in  Figure  6. 
This  analysis  is  one  dimensional  and  assumes  no  axial  (z )  dependence. 

Notice  that  Xe  reaches  a  minimum  at  about  r=1.15  cm  (Xe=8  cm). 
Therefore,  in  the  absence  of  the  magnetic  field,  the  electrons  are  essentially 
collisionless.  It  is  the  presence  of  the  magnetic  field  which  causes  the 
electrons  to  gyrate  about  the  field  lines  in  relatively  small  gyromagnetic 
orbits  (rL  e  ~  5m).  This  promotes  electron-electron  interactions  which  will 
rapidly  provide  a  mechanism  for  heat  transport.  Notice  that  like-particles 
which  collide  (while  executing  gyro-motion)  very  efficiently  transfer 
momentum  and  energy.  The  Larmor  radius  (rLje)  is  given  by: 

_  vthe  _  n  00  Gauss  cm  T* 

©c  eV54  B 
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r  (cm) 


Figure  6:  Electron  Mean  Free  Path  versus  Distance  from  Plasma  Jet 
Axis. 


where  vthe  is  the  electron  thermal  velocity,  (oc  is  the  electron  gyrofrequency. 
The  magnetic  field,  B ,  for  the  purposes  of  these  examinations  is  assumed  to 
be  approximated  by: 

B(r)  =  2309.1  Gauss  cm"2  r2e^,0-2cm  , 
as  is  shown  in  Figure  7a. 

This  profile  ensures  that  the  maximum  magnetic  field  is  50  Gauss  at 
r=0.4  cm  and  0  Gauss  on  the  axis,  r=0  cm.  Using  this,  the  Larmor  radius 
(as  a  function  of  radius)  becomes: 


=  2.38  Gauss  cm  eV 


-vfr((Tmax  ~  ^  +  Tbg? 

2309.1  r2e_r/0'2cm 


This  is  shown  in  Figure  7b.  Notice  that  between  r=0.4  cm  and  r=l.l  cm  the 
Larmor  radius  is  small  enough  so  that  the  electron  will  complete  a  full 
gyration  within  the  size  typical  of  the  plasma  jet.  The  ratio  of  gyrofrequency 
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Figure  7a:  Assumed  magnetic  field  as  a  function  of  range  from  jet  axis. 
This  ensures  B- 0  Gauss  at  r=0  cm  and  B- 50  Gauss  at  r  =0.4  cm. 

to  collision  frequency  for  the  electron  is  given  by: 

(/rp  rp  \  -r2/8^.  .  rp  \3/2 

ooctle  =  0.0152  Gauss-1  eV-3/2  r.u.  - — - — — - - - ^  • 

max  —  nbg  )®  +  nbg  ) 

This  is  shown  in  Figure  8. 

Inside  r=l.l  cm,  the  electrons  are  magnetized  and  execute  gyro-motion 
about  the  magnetic  field  lines.  Although  electron-electron  collisions  in  this 
region  do  not  produce  particle  diffusion,  they  do  efficiently  produce  thermal 
diffusion.  Beyond  r=l.l  cm,  the  electrons  are  unmagnetized  and  move  freely 
with  the  electron  thermal  velocity.  As  shown  in  Figure  6,  the  electrons 
essentially  exist  in  a  collisionless  plasma  (due  to  their  high  temperature) 
and  are  only  constrained  due  to  the  presence  of  a  magnetic  field.  Thermal 
diffusion  of  heat  can  be  described  by  Bohm  diffusion.  An  approximate 
transport  rate  or  diffusion  velocity  can  be  calculated,  and  is  given  by: 
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Figure  7b:  Electron  Larmor  Radius  versus  Distance  from  Plasma  Jet 
Axis. 


Db 


(23) 


where  Dq  is  the  Bohm  diffusion  coefficient  and  5m  is  the  radial  size  of  the 
plasma  jet.  The  Bohm  diffusion  coefficient  is  given  by: 


D=JL?E±-=D 

B  16  eB  1 


Assuming  kTe- 10  eV,  B  =50  Gauss,  and  8™  =0.5  cm;  implies  vb  =4.2x10® 
cm/sec.  Comparing  ug  with  ( v =  1.2x10®  cm/sec)  shows  that  heat 

transport  occurs  faster  than  density  (or  particle)  transport.  This  is 
supported  by  the  program  pl7.c  which  simulates  the  two-dimensional 
diffusion  of  a  plasma  column.  Specific  results  of  pl7.c  will  be  discussed  in 
detail. 
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Figure  8:  Ratio  of  Electron  Gyrofrequency  to  Collision  Frequency 
versus  Distance  from  the  Plasma  Jet  Axis. 

Bohm  diffusion  describes  particle  diffusion  of  a  fully  ionized  plasma  in 
the  presence  of  a  magnetic  field.  Bohm  diffusion  has  historically  provided 
reasonable  agreement  with  experiments.  The  Bohm  diffusion  coefficient  is 
based  on  a  semiempirical  formula.  Notice  that  the  diffusion  coefficient 
depends  on  the  inverse  magnetic  field  and  is  proportional  to  the 
temperature.  This  is  contrary  to  "classical"  cross  field  diffusion  theory, 
where  the  diffusion  coefficient  is  proportional  to  the  inverse  square  of  the 
magnetic  field.  Now  assuming  electron  temperature  diffusion  is  determined 
by  Bohm  diffusion  within  the  electron  fluid,  the  temperature  diffusion 
equation  is: 

-^-  =  v.(z?Bvre),  (25) 

where  the  electron  temperature  is  the  property  being  diffused  across  the 
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magnetic  field.  The  characteristic  temporal  decay  in  Bohm  diffusion  is 
exponential.  The  time  constant  for  this  decay  in  a  cylindrical  plasma  column 
of  radius  (R ),  length  (L),  and  temperature  (T)  can  be  estimated  to  be 
[Ref.4], 


T  _  _  R 

dT/dt  Tr2nRL  pE  2fr  ’ 


where  rr  =  DBdpE/dr,  is  the  energy  flux  and  pE  is  the  energy  density.  Using 
this  in  the  above  equation  gives: 


x 


R‘ 


2 D 


B 


=  XB. 


(26) 


where  xB  is  the  Bohm  time.  For  the  plasma  jet  experiments  examined  here: 
kTe  =  10.0  eV,  B=50  Gauss,  and  R  =0.5  cm;  which  implies  xB  =  100  nsec.  Since 
Texp  (the  characteristic  time  scale  of  the  experiment)  is  the  same  order  of 
magnitude  as  the  Bohm  time,  xB,  detectable  heat  transport  and  radial 
particle  migration  will  occur.  In  agreement  with  experimental  results,  the 
plasma  jet  expands  about  1  cm  in  the  radial  direction  in  the  time  scale  of 
the  experiment.  This  was  supported  by  numerical  simulations  performed  by 
pl7.c,  which  simulated  the  radial  two-dimensional  diffusion  of  a  gaussian 
plasma  column  under  conditions  similar  to  those  encountered 
experimentally.  It  should  be  emphasized  that  the  simulations  shown  in 
Figures  9a  and  9b  show  radial  electron  density  and  temperature  diffusion 
across  the  magnetic  field  assuming  no  axial  flow.  That  is,  the  plasma  is 
assumed  to  be  two  dimensional  (r ,  0)  and  hence  does  not  expand  in  the  axial 
direction.  The  plasma  is  assumed  to  have  azimuthal  symmetry  and  is 
expanding  radially  into  a  magnetic  field  which  is  constant  in  both  space  and 
time.  See  Figures  9a  and  9b.  The  vertical  lines  depicted  in  the  figures 
indicate  that  the  electron  temperature  and  density  are  assumed  to  be 
maintained  constant  by  the  shock  front.  Electron  temperature  and  density 
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are  assumed  to  be  constant  out  to  an  incident  angle  of  60  degrees,  which 
corresponds  to  the  V2  downpoint. 


.  V  U _ I - — 

0  r(cm)  5.0  eT  =  0.002 


Figure  9a:  Numerical  Simulation  of  Thermal  and  Density  Diffusion  us¬ 
ing  pl7.c  at  0.5  nsec.  This  figure  is  enclosed  as  an  example  of  the  ini¬ 
tial  profile  present.  Shock  heating  maintains  electron  temperature  and 
density  constant  out  to  the  vertical  line  (=  0.4  cm). 


A  copy  of  pl7.c  is  enclosed  in  the  appendix.  Further  discussion  about 
simulated  parameters  will  follow  in  subsequent  sections  of  this  thesis.  It  is 
physically  expected  that  the  radial  density  diffusion  will  be  slow  compared 
to  the  time  scales  of  jet  axial  expansion  (10-7  sec).  This  is  because  plasma 
density  diffusion  is  dominated  by  ion  motion  which  occurs  more  slowly,  due 
to  the  slower  ion  sound  velocity.  Figures  9a  and  9b  show  ion  diffusion  and 
thermal  Bohm  diffusion.  Ion  diffusion  will  be  discussed  in  the  next  section. 
The  following  typical  values  were  used  in  the  simulation:  B=50  Gauss, 
kTe=43.6  eV,  and  a  maximum  plasma  density  of  15.0  relative  units.  The 
relative  units  denote  actual  density  measurements  by  electric  double  probes. 
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**  Radial  Cross  Field  Diffusion  ** 

Time=700  nsec  max  n  =  15  r.u. 
'max  : —v  43.597  eV 


8rT/8rn  =  1.3077 
n(4cm)  =  0.10  r.u. 

T(4cm)  =  1.00  eV 
nbg  =  0.10  r.u. 

Tbg  =  1.00  eV 

~x-*-  n(0)  =  15.00  r.u. 

Tft(0)  =  43.60  eV 

8rnQ  =  0.55  cm 

8rT  —  0.55  cm 
o 

B  =  50  gauss 
dt  =  0.001  nsec 

_  en  -  1’621 

r( cm)  5.0  =  1.222 

Figure  9b:  Numerical  Simulation  of  Thermal  and  Density  Diffusion  us¬ 
ing  pl7.c  at  700  nsec.  Shock  heating  maintains  electron  temperature 
and  density  constant  out  to  the  vertical  line  (=  0.4  cm). 

The  relationship  between  the  relative  units  of  density  and  density  in 
particles/cm  3  is  given  by: 

parti^es.  =  3  18xl0nxVmp| 
cm3 

where  Vmp  is  the  measured  potential  drop  in  volts,  as  obtained  by  the 
current  to  an  electric  dipole  density  probe.  This  conversion  factor  can  be 
used  to  convert  to  conventional  density  units;  however,  this  will  prove  to  be 
unnecessary  for  the  work  presented  here.  It  should  be  mentioned  again  that 
axial  expansion  was  intentionally  omitted  (in  pl7.c)  because  the  specific 
effects  due  to  radial  cross  field  diffusion  wanted  to  be  simulated.  This 
objective  was  accomplished. 


T(eV) 


n  (r.u. ) 
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Given  that  the  plasma  jet  is  supersonic  with  Mi  =  11.75,  it  is  now 
possible  to  calculate  (using  equations  (18),  (19),  and  (20))  the  increases  in 
temperature,  pressure,  and  density  across  the  shock.  For  example,  the 
temperature  increases  by  a  factor  of  43.6  across  the  shock  front  assuming 
no  ionization  or  radiation  losses.  The  plasma  is  shock  heated  to 
approximately  43  eV.  Numerical  simulations  (pl7.c)  show  that  when  shock 
heating  occurs,  magnetic  field  reversal  occurs  because  the  radial 
temperature  gradient  relaxes  radially  at  a  faster  rate  than  does  the  density 
gradient.  It  should  be  noted  that  this  phenomenon  also  occurs  in  the  axial 
direction.  Electrons  rapidly  conduct  some  of  their  thermal  energy  outward  in 
front  of  the  shock  and  "pre-heat"  the  incoming  electrons.  This  can  be  seen 
when  comparing  the  shock  velocity,  vsh  =  3xl06  cm/sec,  to  the  Bohm  diffusion 
velocity.  Electron  temperature  propagates  ahead  of  the  shock  and  pre-heats 
the  incoming  electrons.  This  will  be  analyzed  and  discussed  further  in 
Chapter  V  section  F.  Density  diffusion  of  the  plasma  jet  in  the  radial 
direction  across  $  can  be  analyzed  using  classical  particle  diffusion  theory. 
This  is  possible  because  the  ions  are  not  magnetized  (tocxlf  <  1)  and  have  a 
large  mean  free  path  (X*  a>  0.5  cm).  The  ions  essentially  disperse  radially  at 
their  respective  thermal  velocity.  Ion  diffusion  will  be  discussed  further  in 
the  next  section.  If  there  is  no  initial  magnetic  field  present  when  shock 
heating  begins  then  it  is  worth  noting  that  the  electrons  will  not  be 
constrained  by  a  magnetic  field.  Therefore,  the  electrons  will  also  be 
collisionless  due  to  their  high  temperature  and  propagate  away  from  the 
plasma  jet  (within  the  limits  imposed  by  quasineutrality)  and  form  a  bi¬ 
layer.  Chapter  IV  will  show  that  this  freedom  will  generate  even  larger 
shock  generated  fields  (because  is  larger). 

Before  going  on  to  the  next  section  it  should  be  mentioned  that  the 
Tokamak  fusion  reactor  has  produced  similar  "anomalous"  behavior  where 
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the  cross  field  electron  heat  transport  rate  greatly  exceeds  the  particle 
diffusion  rate.  These  large  electron  heat  transport  rates  are  still  under 
study  and  are  not  fully  understood  in  the  fusion  application. 

E.  ION  DIFFUSION 

Particle  diffusion  in  the  radial  direction  is  dictated  by  ion  motion.  In 
the  axial  direction  the  ions  undergo  directed  motion  due  to  the  energy 
received  during  the  laser  blast.  The  mean  free  path  for  the  ions  is  much 
larger  than  the  characteristic  size  of  the  plasma  jet  (X,-  =  28  cm  at  r=0  cm; 
X*  =  2  cm  at  r=l.l  cm).  The  decrease  in  the  mean  free  path  along  r  is  due  to 
the  ion  temperature  gradient.  The  Larmor  radius  for  the  ions  has  a 
minimum  of  about  50  cm  at  approximately  r=0.6  cm.  Therefore,  it  is  evident 
that  the  ions  are  essentially  collisionless  and  travel  straight  paths  since 
their  Larmor  radius  is  so  large  (rLii  »8m).  That  is,  the  ions  are  not 
magnetized  (<ocxlf  <  1).  This  implies  that  the  ions  expand  out  radially  at  the 
ion  sound  velocity  given  by  (assuming  Tt  =  Te ): 

kT  1  **  T  % 

— -  =  9.79xl05  cm  sec-1  eV_V4  -{r  ,  (27) 

mi  J 

where  Te  is  in  eV,  and  |x  is  the  ion  to  proton  mass  ratio,  and  vthi  is  in 
cm/sec.  Using  Te- 43  eV  and  p=28  for  the  shock  heated  plasma  implies 
vthi  =  1.2xl06  cm/sec.  The  shock  velocity  is  3xl06  cm/sec.  Comparing  these 
two  velocities  (vthilvsh  =  0.4)  shows  why  the  plasma  jet  expands  faster  in  the 
axial  direction  than  in  the  radial  direction.  The  jet  experiences  directed 
motion  in  the  axial  direction  while  expanding  at  the  ion  sound  velocity  in 
the  radial  direction.  Given  that  =  1.2xl06  cm/sec,  implies  that  in  200  nsec 
the  plasma  jet  radially  expands  about  0.24  cm  which  is  in  agreement  with 
experiment. 
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Although  electron  heat  transport  occurs  more  rapidly  than  ion  particle 
transport,  electron  particle  transport  occurs  at  the  same  rate  as  ion  particle 
transport.  Due  to  the  coulombic  interaction  between  ions  and  electrons,  any 
charge  separation  results  in  an  electric  field  which  accelerates  electrons 
back  to  the  ions  and  visa  versa.  Therefore,  the  radial  electron  density 
gradients  relax  at  the  ion  sound  velocity.  Classical  particle  diffusion  is 
mathematically  described  by: 

^-=DV2n,  (28) 

where  D  is  the  particle  diffusion  coefficient.  The  relationship  between  the 
particle  diffusion  coefficient  and  other  fundamental  constants  is  given  by 
[Ref.6]: 


where  X*  is  the  mean  free  path  between  collisions  (which  may  also  be 
considered  to  be  the  average  displacement  in  a  collision);  x*  is  the 
characteristic  collision  time,  or  time  between  collisions. 

Using  the  "bi-diffusive”  density-temperature  model  described  in  this 
chapter,  in  conjunction  with  axial  shock  heating,  the  following  chapters  will 
provide  experimental  results  and  numerical  simulations  as  evidence  for  this 
model. 

Shock  generated  magnetic  fields  are  produced  by  the  bi-diffusive 
nature  of  the  plasma,  i.e.  electron  temperature  diffuses  more  quickly  than 
electron  density.  The  shock  front  provides  a  source  of  heating  while  the  bi- 
diffusivity  provides  the  proper  arrangement  of  gradients.  Previous  studies 
completed  on  plasma  jets  provided  several  models  to  specifically  explain 
field  reversal  phenomena.  The  "bi-diffusive”  density-temperature  diffusion 
model  is  the  only  model  consistent  with  source  term  simulations. 
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III.  PREVIOUS  EXPERIMENTAL  WORK 


A,  INTRODUCTION 

Self-generated  magnetic  fields  produced  by  laser  plasmas  were  first 
studied  and  detected  by  Stamper  [Ref.7]  in  1971.  Stamper  used  a 
neodymium  doped  glass  laser  with  an  output  of  60  joules  in  30  nanoseconds 
to  target  a  250  micron  fiber  of  Lucite.  The  Lucite  fiber  was  suspended  in  an 
inert  nitrogen  background.  Magnetic  fields  were  detected  that  expanded 
with  the  same  velocity  as  the  expanding  plasma  jet.  The  fields  were 
essentially  azimuthal  with  respect  to  the  jet.  Magnetic  fields  up  to  1000 
Gauss  were  measured  in  a  nitrogen  background  of  200  mtorr.  Later  in  1971, 
Dean  [Ref.8],  studied  the  interaction  between  two  counterstreaming  plasmas 
produced  by  laser  irradiation  of  fiber  targets  in  a  background  nitrogen  gas. 
Dean  essentially  used  the  same  experimental  apparatus  and  found  the 
following:  (l)a  well  defined  front  exists  at  the  leading  edge  of  the  jet,  (2) 
interactions  which  occur  at  the  front  are  strongly  dependent  on  background 
gas  pressure,  (3)  the  front  thickness  was  approximately  1-2  mm,  (4)  the 
front  was  modeled  as  a  shell  of  increased  density. 

Reversal  of  the  azimuthal  magnetic  field  direction  at  the  leading  tip  of 
the  plasma  jet  at  sufficiently  high  background  pressures  (>200  mtorr)  was 
first  observed  by  McKee  [Ref.l]  at  the  Naval  Postgraduate  School  in  1971. 
McKee  used  a  Korad-k-1500  Q  switched  neodymium  doped  glass  laser  which 
delivers  approximately  6.6  Joules  in  a  pulse  of  22  nanosecond  duration.  This 
equates  to  approximately  300  megawatts  of  power.  These  intense  pulses  of 
laser  light  then  irradiate  aluminum  targets  which  have  been  placed  in  an 


31 


evacuated  chamber  typically  between  5x1  CT5  and  0.7  torr  of  inert  gas 
pressure.  The  laser  strikes  the  target  at  30  degrees  off  normal  to  prevent 
laser-target  plasma  interference.  Since  that  discovery,  several  graduate 
students  at  the  Naval  Postgraduate  School  have  attempted  to  explain  this 
reversal  phenomenon.  This  is  the  first  thesis  known  to  the  author  which 
proposes  that  this  reversal  is  the  fingerprint  of  shock  generated  magnetic 
fields.  Previous  experimental  works  will  be  discussed  in  this  chapter. 

B.  MCKEE’S  "SHELL"  MODEL 

McKee  [Ref.l]  later  adopted  the  "shell"  model  in  his  Ph.D.  thesis  to 
account  for  the  phenomenon  of  magnetic  field  reversal  at  the  jet  front  at 
later  times  (300  nsec)  in  background  pressures  greater  than  200  mtorr.  In 
the  "shell"  model  the  expanding  plasma  jet  is  thought  to  cause  a  "pile  up"  of 
background  gas  at  the  leading  edge  of  the  jet.  This  "pile  up"  forms  a  region 
of  high  density  directly  in  front  of  the  plasma  jet.  The  expanding  laser 
plasma  acts  to  compress  and  heat  the  background  in  a  "shell"  which 
surrounds  the  expanding  front  of  the  jet.  See  Figure  10. 

This  "shell"  travels  with  the  jet  and  also  causes  additional 
temperature  and  density  gradients  to  be  formed.  In  this  way,  the 
background  gas  interacts  with  the  expanding  laser  plasma.  McKee 
additionally  attributed  field  reversal  at  background  pressures  higher  than 
200  mtorr  to  these  secondary  gradients  formed  by  the  "shell".  That  is,  as  the 
expanding  jet  formed  a  frontal  "shell",  the  region  between  the  shell  and  the 
laser  plasma  would  develop  temperature  and  density  gradients  in  the  +£ 
direction  (see  Figure  10).  The  source  term  in  equation  (10)  shows  that  if  the 
direction  of  the  axial  gradients  reverses,  then  the  cross  product  changes  sign 
and  magnetic  field  reversal  begins.  This  is  seen  more  explicitly  if  the  source 
term  is  written  as: 
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Figure  10:  Shell  Model  proposed  by  McKee.  [Ref.2] 
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where  now  the  axial  gradients  are  in  the  +£  direction.  In  the  shell  model, 
the  shell  never  actually  transverses  the  plasma  jet  front,  and  therefore 
differs  greatly  from  the  model  proposed  in  this  thesis.  McKee’s  work 
provided  extensive  spatial  and  temporal  two-dimensional  mappings  of  the 
magnetic  fields  at  three  principle  background  pressures:  0.1  mtorr,  5  mtorr, 
and  250  mtorr.  McKee’s  data  provided  very  little  contour  data  concerning 
reversal  phenomena  since  this  phenomenon  is  not  discernible  until  about 
200  mtorr.  However,  some  reversal  effects  are  seen  at  250  mtorr.  McKee 
provided  the  following  conclusions:  (1)  the  self-generated  magnetic  fields 
were  indeed  generated  by  conduction  currents  producted  by  non-parallel 
gradients,  (2)  fields  were  initially  clockwise  as  viewed  in  Figure  1,  then  as 
field  reversal  occurs  the  fields  are  counterclockwise  (Figure  3),  (3)  the  fields 
are  azimuthally  symmetric,  and  (4)  the  magnitudes  of  the  generated 


33 


magnetic  fields  are  amplified  by  the  existence  of  a  background  gas. 
According  to  McKee’s  conclusions,  the  initial  field  direction  implies  that: 


dTe  dne 

dTe  dne 

dr  dz 

> 

dz  dr 

where  when  shock  heating  and  hence  field  reversal  occurs: 


dTe  dne 

dTe  dne 

dz  dr 

> 

dr  dz 

C.  BIRD’S  WORK 

Bird  [Ref.2]  continued  McKee’s  investigation  of  self-generated  magnetic 
fields  produced  by  laser  plasmas.  McKee’s  study  was  mainly  limited  to 
magnetic  contour  mapping  and  characteristics  of  the  fields  in  various 
background  pressures  of  nitrogen  up  to  250  mtorr.  Bird  completed  a 
thorough  plasma  density  mapping  using  electric  double  probes.  The 
additional  data  provided  by  Bird  made  it  possible  to  associate  particular 
magnetic  phenomena  (such  as  field  reversal)  with  specific  density  structures 
of  the  laser  plasma.  Bird  also  analyzed  the  behavior  of  plasma  jets  in 
various  background  gases  (such  as  H2,  He,N2,  and  Ar)  up  to  700  mtorr. 
This  provided  good  field  reversal  data.  Bird  also  conducted  some  relatively 
high  pressure  experiments  at  5  torr  H2.  The  results  are  shown  in  Figure  11 
which  shows  the  measured  density  profiles  n*  and  the  magnetic  field  B  as 
functions  of  time  at  locations  r=0.4  cm  from  the  axis,  0  =  0°,  and  z  =0.4  cm, 
0.8  cm,  1.0  cm,  and  2.0  cm  from  the  target  surface.  Field  reversal  is 
dominant  at  the  jet  front.  In  fact,  the  reverse  field  magnitude  grows  very 
rapidly  until  it  is  the  dominant  field  present.  This  high  pressure  example 
shows  that  the  shock  mechanism  produces  significant  reversed  fields. 
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Bird  was  also  able  to  draw  definite  conclusions  about  the  spatial 
relationships  between  the  generated  magnetic  fields  and  measured  plasma 
jet  density  distribution.  Some  of  Bird’s  major  conclusions  follow.  First,  the 
"shell"  model  proposed  by  McKee  to  explain  the  background  pressure 
dependence  of  the  magnetic  field  was  flawed.  Density  mappings  completed 
by  Bird  were  unable  to  detect  a  density  "shell".  Work  performed  by  Dean 
[Ref.8]  showed,  using  three  independent  techniques,  that  a  "shell"  type 
structure  did  exist;  therefore  Bird  concluded  that  an  analogous  structure 
must  actually  exist  at  the  front  of  the  expanding  plasma  jet.  The  density 
"shell"  found  by  Dean,  however,  is  different  than  the  "shell"  conceived  by 
McKee  in  Figure  10.  It  is  contended  in  this  thesis  that  the  "shell"  structure 
found  by  Dean  is  instead  the  shock  front  at  the  leading  edge  of  the  plasma 
jet.  Second,  Bird  proposed  a  model  which  attributes  the  increase  in  the 
magnetic  field  with  background  pressure  to  interactions  between  the 
photoionized  background  plasma  and  laser  plasma  during  laser  irradiation. 
Third,  because  Bird  performed  his  investigation  at  background  pressures 
above  200  mtorr,  he  observed  extensive  reverse  field  phenomena.  Bird 
observed  that  field  reversal  appeared  to  be  due  to  a  "pile  up  or  (re- 
thermaHzation)"  of  the  laser  plasma  as  it  propagated  through  the 
background  plasma.  This  "pile  up"  produces  axial  electron  temperature 
gradients  at  the  laser  plasma  front  and  reverses  the  generated  magnetic 
field  in  accordance  with  the  source  term  in  equation  (10).  Bird,  however, 
explains  the  heating  mechanism  at  the  "pile  up"  to  be  convective  in  nature 
(again,  this  thesis  contends  that  the  heating  mechanism  is  instead 
completely  attributable  to  shock  heating).  Fourth,  reverse  fields  are 
generated  long  after  the  cessation  of  the  laser  pulse  (1  (xsec).  Fifth,  the  onset 
of  field  reversal  is  delayed  at  lower  background  pressures  because  "pile  up" 
of  the  background  plasma  requires  more  time.  This  effect  can  be  explained 
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using  the  shock  model.  Since  shock  heating  is  the  mechanism  for  field 
reversal,  at  lower  background  pressures  weaker  shock  fronts  are  formed.  A 
weaker  shock  front  provides  less  shock  heating  and  a  smaller  discontinuity 
in  density  and  temperature,  hence  density  and  temperature  gradients  are 
smaller.  Reduced  gradients  lengthen  the  onset  of  field  reversal  and  also 
decreases  the  overall  magnitude  of  the  fields  produced.  Bird’s  density  data 
was  used  to  produce  a  mathematical  model  of  the  density  mappings.  These 
mathematically  modeled  density  profiles  were  then  used  in  a  computer 
simulation  to  predict  generated  magnetic  fields  produced  by  shock  heating. 

D.  BROOKS’  WORK 

Brooks  [Ref.9]  in  late  1973  conducted  laser  plasma  experiments, 
similar  to  McKee  and  Bird,  in  a  vacuum  of  2.5xl0'6  torr  (air).  Brooks 
provided  detailed  analysis  of  plasma  jet  free  expansion  and  target  damage 
due  to  irradiation.  Two  dimensional  contour  mappings  of  plasma  density 
were  created  which  confirmed  the  fact  that  density  gradients  begin  to  decay 
from  the  onset  of  jet  expansion.  Relaxation  times  for  radial  gradients  were 
actually  much  longer  than  relaxation  times  associated  with  axial  gradients. 
Brooks  explained  this  to  be  due  to  the  fact  that  a  radial  "pinching"  (JxZ?) 
force  prevented  expansion  in  the  radial  direction.  A  maximum  field  of  200 
Gauss  corresponded  to  a  current  density  of  about  400  amps  /cm2.  This 
indicated  the  presence  of  a  JxZ?  force  density  of  approximately  7xl04  N/m3. 
The  ratio  of  the  JxB*  force  to  the  radial  pressure  gradient  force  density  is 
less  than  one.  The  radial  pressure  gradient  force  is  given  by: 

=  6.6x10®  N/m3 

where  /  is  the  characteristic  length  of  the  density  gradient  (~  0.5  cm ), 
n  =  4.8xl013  cm-3,  and  assuming  a  conservative  kT  =43  eV.  That  is,  the 
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"pinching"  alone  does  not  prevent  expansion  in  the  radial  direction.  The 
relatively  slow  ion  thermal  velocity  provides  "containment"  in  the  radial 
direction  (see  Figure  12).  Notice  the  free  expansion  in  the  axial  direction. 

A  plot  of  « Txfi  force  density  measured  by  McKee  at  120  nsec  in  a 
background  of  N2  gas  (250  mtorr)  is  shown  in  Figure  13.  Arrow  length  is 
proportional  to  force  density.  The  largest  force  density  at  this  time  was 
1.4xl04  N/m3.  Brooks  concluded  that  the  JxZ?  force  provided  at  least  "partial 
containment"  of  the  laser  plasma;  however,  he  never  hypothesized  what 
other  confinement  mechanisms  may  be  responsible.  An  additional 
contributor  to  the  slow  radial  expansion  is  the  directed  motion  of  the  ions  in 
the  axial  direction  due  to  the  laser  blast  effects. 

E.  BARRIER  SHOCK 

Schwirzke  [Ref.  10, 11]  published  results  in  the  summer  of  1973  which 
showed  that  if  a  glass  plate  is  placed  in  the  path  of  the  expanding  jet,  then 
the  rapid  deceleration  of  the  jet  at  the  plate  causes  a  shock  to  form.  This 
shock  also  causes  field  reversal,  but  at  pressures  far  below  the  pressures 
required  to  achieve  shock  heating  without  an  impinging  glass  plate  barrier. 
Figure  14  shows  the  results  of  an  experiment  performed  at  5  mtorr  N2  with 
and  without  an  obstructing  glass  plate.  The  glass  plate  was  placed  at  z  = 
1.15  cm,  and  data  was  recorded  along  r  =  0.3  cm.  Data  for  the  no  plate  case 
was  recorded  at  z  -  1.0  cm,  r  =  0.3  cm.  Notice  that  up  to  100  nsec  the 
generated  magnetic  field  for  the  two  cases  is  the  same.  After  100  nsec 
however,  the  reversed  field  component  is  clearly  distinguishable. 
Propagating  1.15  cm  in  100  nsec  corresponds  to  a  velocity  of  1.15xl07  cm/sec, 
which  agrees  well  with  the  measured  jet  velocity  at  this  background 
pressure.  After  700  nsec  the  two  cases  again  approach  the  same  field  value. 
This  shows  that  field  reversal  can  be  achieved  at  low  background  pressures 
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Figure  12:  Density  Contours  of  the  Free  Expansion  of  the  Plasma  Jet 
at  Time  200  nsec.[Kef.91 


Figure  13 :  Jxlf  Force  Exerted  of  a  Plasma  Jet.  force  density  at 
120  nsec  for  250  mtorr  of  N2  background  gas.  The  magnitude  of  the 
force  at  this  time  is  1.4xl04  N/m3.  [Ref.l] 

if  the  plasma  jet  is  rapidly  decelerated  by  a  barrier.  The  rapid  deceleration 
causes  shock  formation  and  hence  generation  of  a  shock  magnetic  field 
whose  azimuthal  component  is  reversed  in  direction  near  the  shock  location. 
The  bottom  curve  in  Figure  14  shows  the  magnitude  of  the  reversed  (shock 
generated)  field  alone.  This  was  derived  by  subtracting  the  no  plate  curve 
from  the  plate  curve.  Notice  the  maximum  reversed  (shock  generated)  field 
approaches  200  Gauss  at  about  340  nsec. 
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IV.  NUMERICAL  SIMULATION 


A  PL2.C 


As  mentioned  earlier,  a  numerical  computer  simulation  was  written  to 
aid  in  understanding  which  arrangement  of  density  and  temperature 
gradients  yielded  the  observed  phenomena.  In  this  way  it  became  possible  to 
indirectly  obtain  estimates  of  the  required  electron  temperature  gradients 
which  were,  before  this  thesis,  unknown  and  unmeasured.  Bird’s  work 
[Ref.2]  provided  excellent  data  to  create  a  mathematical  description  of  the 
density  gradients  present  at  high  background  pressures.  Once  mathematical 
expressions  for  density  and  temperature  gradients  were  obtained,  a 
simulation  calculated  VTexVne,  the  field  production  rate.  Section  C  discusses 
the  actual  structure  of  the  program.  This  section  will  motivate  the  purpose 
of  the  program.  A  surprising  result  showed  that  magnetic  field  reversal 
occurred  under  only  one  combination  of  the  temperature  and  density 
gradients.  Let  denote  the  characteristic  axial  gradient  width  of  the 
temperature,  and  5m  denote  the  characteristic  radial  gradient  width  of  the 
tensity.  Reverse  (shock  generated)  magnetic  field  production  at  the  front  of 
plasma  jet  requires  that: 


lV2TellVrne 


ivrreiiv2ne 


>  1. 


Approximating  the  gradients  by: 


n 


T. 


VzTe  ~  - —  ,  V2 ne  ~  - —  ,  WrTe  -  -z —  ,  Vrne  = 


JzT 


drT 


ne 

5 m’ 


42 


gives: 

MW.  i 
(5,r)(5m)  • 

Shock  compression  in  the  axial  direction  implies  that  the  shock  width  is: 

&zT  ~  $zn  while  52t  «  5rr* 

Rapid  electron  heat  conduction  in  the  radial  direction  ensures  that: 

Srr  >  Sm. 

Therefore,  the  ratio  of  the  characteristic  widths  given  above  is  indeed 
greater  than  one,  and  reversed  fields  are  produced.  Reverse  fields  will 
hereafter  be  referred  to  as  shock  generated  fields.  As  mentioned  earlier,  in 
actuality  electron  pre-heating  occurs  in  the  axial  direction.  However,  the 
pre-heating  effect  does  not  lengthen  Szr  to  the  extent  that  it  causes  the 
above  ratio  to  be  less  than  one.  The  radial  temperature  diffusion  is  the 
dominant  mechanism. 

Although  density  gradients  were  relatively  well  known  by  measured 
results,  temperature  gradients  were  adjusted  to  yield  the  observed 
phenomena.  It  quickly  became  apparent  that  maximum  electron 
temperatures  on  the  order  of  40  eV  were  sufficient  to  produce  the  observed 
magnetic  fields.  It  should  be  noted  that  the  simulation  does  not  calculate 
the  magnetic  fields  generated  very  early  on  (<30  nsec)  by  the  incident  laser 
radiation.  At  these  early  times,  magnetic  field  production  is  dominated  by 
radial  temperature  gradients  created  by  the  laser  and  axial  density 
gradients  created  by  the  expanding,  ionized,  dense  target  material.  As 
discussed  earlier,  these  fields  diffuse  slowly  and  tend  to  be  carried  out  with 
the  plasma  jet.  After  cessation  of  the  laser,  the  "frozen  in"  magnetic  field  is 
modified  by  the  shock-generated  magnetic  fields  created  by  the  gradients 
present.  The  numerical  simulation,  pl2.c,  simulates  the  post-laser  self- 
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generated  magnetic  fields  only.  pl2.c  is  a  computer  code  written  in  the  C 
programming  language.  A  copy  of  pl2.c  is  attached  as  an  appendix  at  the 
end  of  this  thesis.  The  density  profiles  modeled  by  pl2.c  were  obtained  from 
experimental  data  taken  by  Bird.  To  ensure  shock  formation  and  hence  field 
reversal,  pl2.c  uses  data  obtained  from  measurements  taken  at  a 
background  pressure  of  700  mtorr.  Other  pressures  could  have  easily  been 
used,  however,  700  mtorr  provides  the  most  interesting  example.  Figures 
15a  through  15h  show  Bird’s  data  for  this  case.  All  these  measurements 
were  made  at  r=0.4  cm  and  0  =  0°.  This  implies  that  the  data  does  not 
accurately  represent  the  plasma  density  present  for  times  less  than  200 
nsec  because  the  plasma  front  does  not  arrive  until  that  time.  This  can  be 
seen  when  comparing  Figures  15b  and  15c.  The  shock  generated  magnetic 
field  shows  up  at  300  nsec  at  the  steepest  part  of  the  axial  density  profile. 
For  a  given  set  of  gradients,  the  magnetic  field  increases  with  time: 

f*  dt  =  B(*2)  -  ( VTexVne)dt  . 

1  ut  6/lg 

B.  DENSITY/TEMPERATURE  PROFILES 

The  mathematical  density  profiles  used  in  pl2.c  assume  an  initial 
gaussian  radial  distribution  and  an  exponential  axial  distribution. 
Representation  of  the  radial  plasma  density  profile  by  a  gaussian  provided 
an  easy  and  fairly  accurate  model.  Axial  density  behavior  was  modeled  by 
an  exponential  which  provided  a  rather  accurate  representation  of  the  shock 
front.  Within  0.3  cm  (the  approximate  shock  thickness)  the  density  rises 
from  background  density  values  (nbg)  of  about  0.2  relative  units  to  the 
measured  results  of  about  15  relative  units.  A  typical  density  profile  fit  to 
experimental  data  has  the  form: 

n(rj,t)  =  (nmex(t)-nbgXl-e(* -*o(tMm)e~r  nrH+nbg  ,  (30) 
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Figure  15a:  B-Field  and  Density  Data  at  20  nsec  in  700  mtorr  N2  along 
a  line  r=0.4  cm  and  8  =  0°.  Solid  line  indicates  plasma  density.  Dashed 
line  indicates  magnetic  field  magnitude.  [Ref.2] 

where  n  max(t )  is  a  time  dependent  function  that  describes  the  maximum 
plasma  density.  z0(t)  is  also  a  time  dependent  function  which  describes  the 
location  of  the  plasma  jet  front.  In  pl2.c  it  is  assumed  that  z0(t )  is  a  linear 
function  of  time.  That  is,  no  accelerations  of  the  jet  front  occur.  The  factor 
in  the  exponent  ensures  that  the  shock  front  has  a  thickness  of  about  3 
mm.  Comparison  with  Bird’s  data  in  Figures  15a  through  15h  shows  that 
the  front  thickness  is  such  that  the  density  increases  to  90%  of  its 
maximum  value  within  5  mm.  This  observation  was  then  used  to  determine 
the  exponential  axial  constant  of  5^  =  0.25  cm.  The  gaussian  factor  describes 
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Figure  15b:  B-Field  and  Density  Data  at  60  nsec  in  700  mtorr  N2  along 
a  line  r  =0.4  cm  and  0  =  0°.  Solid  line  indicates  plasma  density.  Dashed 
line  indicates  magnetic  field  magnitude.  [Ref.2] 

the  radial  behavior  of  the  plasma.  The  factor  in  the  exponent  indicates 
that  at  r  =  5m  cm,  the  density  has  fallen  off  by  a  factor  of  e-1  or  63%.  It  was 
assumed  that  the  characteristic  radial  width  of  the  density  gradient  8m 
remained  constant.  This  assumption  is  supported  by  Figures  9a  and  9b.  It  is 
evident  from  these  figures  that  the  radial  growth  is  minimal.  The  second 
term  of  equation  (30)  (n^=  0.2  r.u.)  represents  the  background  plasma 
density.  A  value  of  0.2  relative  units  corresponds  to  about  6.4xl0n  cm  3. 
Figure  16  shows  the  experimentally  measured  relationship  between  the 
maximum  jet  density  versus  time.  Notice  that  nmax  versus  time  is  fairly 
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Figure  15c:  B-Field  and  Density  Data  at  100  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  0  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

independent  of  background  pressure.  It  is  again  important  to  notice  that 
this  data  was  taken  at  r=0.4  cm,  0  =  0°.  The  plasma  jet  does  not  arrive  at 
the  measuring  probe  until  some  time  between  200  nsec  and  300  nsec.  This 
is  easily  seen  in  Figure  16.  Data  for  background  pressures  of  0.5  mtorr  and 
5  mtorr  were  only  recorded  until  300  nsec.  This  is  because  all  relevant 
magnetic  field  activity  had  already  decayed  away  at  these  lower  pressures. 

The  functional  dependence  of  nmax  at  r=0  was  extrapolated.  The 
mathematical  form  of  nmax(t)  ensures  that  nmax( 300)  =  15  and  nmax( 700)  =  5 
which  matches  the  extrapolated  experimental  data.  z0(t)  is  assumed  to  be  a 
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Figure  15d:  B-Field  and  Density  Data  at  200  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  9  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

linear  function  of  time.  This  is  supported  by  Brooks’  data  shown  in  Figure 
17.  Plasma  front  velocity  is  a  constant.  At  700  mtorr  N  2  background 
pressure,  the  jet  velocity  was  about  3.1xl06  cm/sec,  hence  z0{t )  is  given  by: 

z0(t )  =  3.1xl0'3  cm/nsec  t  , 

where  t  is  in  nsec  and  z0  is  in  cm.  In  1000  nsec  the  jet  will  propagate  about 
3.1  cm.  The  graphical  representation  of  n(r^,t)  is  shown  in  Figure  18a. 

The  electron  temperature  is  modeled  by  an  equation  very  similar  to 
equation  (30).  No  experimental  results  exist  which  provide  data  on 
temperature  distributions  within  the  vicinity  of  the  plasma  jet,  therefore  it 
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Figure  15e:  B-Field  and  Density  Data  at  300  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  0  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

was  assumed  that  the  radial  distribution  was  also  a  gaussian  and  that  the 
axial  distribution  was  exponential.  Curvature  of  the  plasma  jet  front  implies 
that  the  centerline  temperature  of  the  jet  is  greater  than  the  "off  axis" 
temperature.  This  is  due  to  the  fact  that  shock  heating  occurs  only  to  the 
fluid  elements  whose  velocity  components  are  normal  to  the  shock  front.  An 
exponential  axial  distribution  implies  that  the  jet  front  not  only  provides  a 
rapid  rise  in  plasma  density,  but  also  in  plasma  temperature.  Shock  heating 
at  the  jet  front  generates  this  rise  in  electron  temperature.  A  typical 
mathematical  temperature  profile  (fit  to  numerical  results)  is  given  by: 
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Figure  15f:  B-Field  and  Density  Data  at  400  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  0  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

T(r^,t)  =  (Tmax-TbgXl-eiz~ZoUW^T)e~r2l5?r{t)  +  Tbg  ,  (31) 

where  T(r#,t)  and  Tmax  are  in  eV,  5zT  =  1/3.99  cm,  SrT  is  the  characteristic 
radial  width  of  the  temperature  gradient,  and  Tbg  is  the  electron 
background  temperature.  The  exponential  factor  (§rr)  in  the  above 
expression  for  the  temperature  profile  is  not  exactly  equal  to  0.25  cm  (as  in 
the  case  of  the  density  profile)  due  to  the  effects  of  axial  electron  heat 
diffusion  "pre-heating"  incoming  electrons.  Since  electron  heat  conduction 
occurs  very  rapidly,  hrT{t)  will  not  be  a  constant  function  of  time.  The 
mathematical  form  of  T(r^,t)  is  shown  in  Figure  18b.  Using  the  forms  of 
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Figure  15g:  B-Field  and  Density  Data  at  500  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  0  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

T(r^,t)  and  n(rj,t)  given  above,  pl2.c  solves  the  time  integrated  source 
equation  (10).  Magnetic  field  diffusion  and  convection  are  not  calculated  by 
pl2.c. 


C.  FINITE  DIFFERENCE  EQUATIONS 

The  pl2.c  program  simulates  the  production  of  magnetic  flux  due  to  the 
source  term  in  equation  (10).  The  source  term  is: 


dff 

dt 


=  -VTe  x  — Vne  = 
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Figure  15h:  B-Field  and  Density  Data  at  700  nsec  in  700  mtorr  N2 
along  a  line  r=0.4  cm  and  6  =  0°.  Solid  line  indicates  plasma  density. 
Dashed  line  indicates  magnetic  field  magnitude.  [Ref.2] 

pl2.c  performs  a  point  by  point  time  integration  to  determine  the  generated 
magnetic  field  at  every  spatial  mesh  point.  Given  that  the  plasma  jet 
typically  has  dimensions  of  1.5  cm  in  diameter  and  about  4  cm  in  length 
prior  to  dissipating,  this  region  of  space  is  divided  into  a  250  by  250  two- 
dimensional  point  grid.  Each  point  on  the  grid  is  assigned  a  coordinate 
(z,r),  a  temperature  T(z,r),  a  density  n(z,r),  and  a  magnetic  field  B(zf). 
The  radial  coordinate  (r  )  varies  from  -1.5  cm  to  1.5  cm.  In  grid  coordinates, 
the  radial  coordinate  is  denoted  by  an  integer  j  which  varies  from  0  to  250. 
The  jet  expands  along  the  j= 125  axis  which  corresponds  to  the  r=0  cm  axis. 
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Figure  16:  nmax  versus  Time  (Experimental)  at  r=0.4  cm,  0  =  0°.  z 
varies  depending  on  where  the  maximum  density  occurs.[Ref.2] 


The  axial  coordinate  (z)  varies  from  0  to  twice  the  jet  front  expansion 
distance.  That  is,  if  the  jet  expands  to  1.5  cm,  then  the  grid  varies  from  0  to 
3.0  cm.  The  axial  grid  coordinate  is  denoted  by  the  integer  i,  and  (as  j) 
varies  from  0  to  250.  The  grid  coordinates  then  become  (i ,  j )  and  the 
temperature  and  density  profiles  become  point  defined  at  each  of  these  L 
coordinates.  Grid  profiles  are  denoted  by  Tij,  riij,  and  Btj.  In  order  to 
numerically  calculate  magnetic  field  growth  at  a  point,  equation  (32)  must 
be  converted  to  a  finite  difference  equation  [Ref.  12].  That  is,  it  must  be  cast 
into  a  form  which  can  be  utilized  by  the  computer.  The  partial  derivatives  in 
equation  (32)  at  point  (i,j)  can  be  calculated  in  terms  of  the  neighboring 
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Figure  17:  Plasma  Jet  Velocity.  Position  of  the  front  of  the  plasma  jet 
versus  time  for  various  background  gas  pressures.  Data  was  taken 
along  r  -  0  cm  with  Langmuir  double  probe.  [Ref.9] 


grid  points  (i+l,  j),  (i,j+ 1),  (i-1,  j),  and  (i,j- 1).  The  partial  derivative  of  the 
electron  temperature  with  respect  to  the  axial  coordinate  0 Te/dz),  for 
example,  is  written  in  terms  of  the  finite  difference  equation  as: 


dTdz  = 


Tj+ij 
2 dz 


9 


where  dTdz  denotes  dTe/dz,  and  dz  is  the  finite  distance  between  the  grid 
points  (i+l,j)  and  (i-l,j).  In  pl2.c:  dz=t  1/250  cm,  where  1 1  is  the 
distance  in  cm  spanned  by  the  250  point  axial  grid.  Since  this  distance 
depends  on  the  integration  time  of  the  simulation,  1 1  =  6.2xl0-3  tmax.  tmax  is 
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Figure  18a:  Mathematical  Density  Profile,  i.e.  n(z,r,t).  See  equation  30. 


the  total  time  (in  nsec)  of  the  simulation.  Notice  that  1 1  corresponds  to  twice 
the  jet  expansion  distance.  The  above  finite  difference  equation  essentially 
calculates  the  slope  of  the  temperature  grid  in  the  z  direction.  Likewise: 


dndz  = 

dTdr  = 


ni+lJ  ni-lJ 
2dz 

Tjj+i  ~  Tjj- 1 
2 dr 


and: 


dndr  = 


j+i  j- 1 

2 dr 


y 


where  dr  =  3/250  cm,  since  the  total  distance  spanned  in  the  radial  direction 
is  always  3  cm. 
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Figure  18b:  Mathematical  Temperature  Profile,  i.e.  T(z,r,t).  See  equa¬ 
tion  31. 


For  each  time  step  ( dt )  the  density  and  temperature  gradients  are 
calculated  for  each  grid  point  with  the  equations  given  above.  Once  the 
partials  are  calculated,  the  magnetic  field  must  be  integrated  at  each  grid 
point.  That  is,  if  equation  (32)  is  written  in  finite  difference  form: 


&ij  new  —  j  0id  +  — —  (dTdz  )(dndr )  -  (dTdr  Xdndz )]  dt  ,  (33) 

enij  L  J 


then  the  magnetic  field  at  (i,j)  is  incremented  by  the  second  term  in  the 
above  equation.  BiJ0id  represents  the  "old"  value  of  the  magnetic  field  at 
grid  point  (i ,j)  at  time  t  -  dt.  The  above  expression  is  nested  in  a  time  loop, 
so  that  the  last  Bijnew  becomes  the  present  Bij0id.  The  magnetic  field 
throughout  the  grid  is  zero  at  t=  0  nsec,  and  grows  as  the  second  term 
increments  the  previous  time  step  value  of  Btj  by  an  expression  related  to 
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the  density  and  temperature  gradients.  Since  the  total  integration  time  is 
typically  on  the  order  of  about  700  nsec,  dt  is  chosen  such  that  at  least  1000 
time  steps  are  completed  in  order  to  calculate  the  magnetic  field.  For  the 
simulations  run  by  pl2.c,  dt  =  0.02  nsec.  The  experimental  values  of  B ,  kT , 
time,  and  length  are  in  units  of  Gauss,  eV,  nsec,  and  cm,  respectively.  In 
order  for  the  calculated  magnetic  field  to  be  in  units  of  Gauss,  temperature 
(kT)  to  be  in  units  of  eV,  time  to  be  in  units  of  nsec,  a  conversion  factor 
must  be  introduced  to  the  second  term  to  ensure  the  proper  outcome  of 
units.  The  density  is  expressed  in  relative  units  since  the  density  units 
cancel  out  in  the  second  term  of  equation  (33).  The  conversion  is: 

(1.6xl0~19  -4r)  GO-9  -^-)  (104  -£-) 

_ eV  nsec _ T_  _  iq-1 

(1.6xl0~19  C )  (10-2  — )2 
cm 

Now  equation  (33)  is  written: 

BtJ  _  =  W  ^  [« *  Xdndr)  -  ' 

This  is  the  form  of  the  source  equation  used  in  pl2.c.  In  pl2.c,  the  point 
functions  nfj,  Ttj,  and  Btj  are  represented  as  two  dimensional  arrays 
denoted  by  n[nz][nr],  T[nz][nr],  and  B[nz][nr],  respectively,  nz  and  nr  are 
also  integers  (0<nr<250,  0<nz<250)  which  represent  the  grid  coordinates  of  a 
point  in  space. 

In  order  to  allow  the  jet  to  propagate  away  from  the  target  surface, 
pl2.c  propagates  a  "nulling"  plane  starting  at  the  jet  front  (at  t=0  nsec)  and 
ending  at  the  target  surface  at  t=tmax.  That  is,  the  plasma  jet  is  stationary 
and  the  "nulling"  plane,  which  can  be  thought  of  as  the  target  surface, 
propagates  away  from  the  plasma  front.  The  jet  frame  of  reference  provides 
the  same  results  as  the  laboratory  frame.  The  jet  frame,  however,  is  a  more 
convenient  method  to  propagate  the  jet.  This  is  simply  an  artifact  of  the 
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simulation  and  introduces  no  additional  error  into  the  analysis. 

D.  PL7.C 

In  order  to  achieve  an  understanding  of  the  radial  diffusion  of  the 
plasma  density  and  temperature,  pl7.c  was  written.  pl7.c  numerically  solves 
the  coupled  partial  differential  equations  describing  temperature  and 
density  diffusion.  As  mentioned  earlier,  the  radial  density  diffusion  is 
modeled  by  classical  particle  diffusion.  The  radial  temperature  diffusion  is 
modeled  by  a  semi-empirically  derived  Bohm  diffusion  coefficient  which 
describes  diffusion  across  a  magnetic  field  in  a  plasma  which  without  a 
magnetic  field  would  be  collisionless.  The  Bohm  diffusion  coefficient  (DB)  is 
given  by  equation  (24)  and  is  itself  a  function  of  electron  temperature  and 
the  magnetic  field  present.  Equation  (25)  represents  the  diffusion  equation 
incorporating  the  Bohm  coefficient  as  an  upper  bound.  The  classical  particle 
diffusion  coefficient,  D,  is  a  function  of  the  ion  collision  frequency  and  ion 
mean  free  path  (given  by  equation  (29)).  That  is,  it  is  also  a  function  of 
temperature  and  density.  Therefore,  pl7.c  solves  the  coupled,  partial 
differential,  diffusion  equations  (25)  and  (28). 

pl7.c  assumes  axial  and  azimuthal  symmetry  and  therefore  reduces  to 
a  one  dimensional  simulation.  Equations  (25)  and  (28)  are  converted  to 
finite  difference  form,  yielding: 

rr  _  rr  .  Mrj+1(TJ+1  Tj+ 2  -  Tj+i  Tj )  +  rj_i  ( Tj _i  Tj_2  -  T/_i  Tj ))  dt 

4.0O  dr2  ’ 


>,j  ~  nold,j  + 


C(rj+ 1  (nj+2  -  rij)  +  rj_x  (n;_2  -  nj))  dt 


4.0/v  dr2 


respectively.  The  integration  constants  A  and  C  are  given  by, 
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A  = 


6.25X10"3 

1 

1 

104g 

1.6x10_19J 

' 

10-9sec 

100cm 

B 

16B 

1.6xl0-19  C 

IT 

leV 

1  nsec 
••  ^ 

lm 

*  d 

and 


4.253xl0"6  T?  5 


where  B  is  the  magnetic  field  in  Gauss  and  rij  is  the  particle  density  in  r.u. 
The  above  finite  difference  equations  are  then  integrated  by  the  Euler- 
Cromer  method.  Values  of  the  other  constants  are:  dt  =0.001  nsec  and 
dr  =3.0/125.0  cm.  The  initial  temperature  and  density  distributions  are  again 
assumed  to  be  gaussian  with  their  respective  temperature  and  density 
widths,  5 rT  and  8m ,  initially  assumed  to  be  equal.  Initial  and  final  total 
mass  and  energy  are  also  calculated.  In  the  early  stages  of  programming, 
this  was  done  to  ensure  simulation  accuracy;  however,  at  later  stages  this 
calculation  was  used  to  determine  how  much  energy  and  mass  was  being 
introduced  due  to  shock  heating.  Shock  heating  ensures  that  the  internal  jet 
temperature  and  density  remain  at  shock  determined  values.  Therefore, 
pl7.c  maintains  the  temperature  and  density  constant  for  the  duration  of 
the  simulation  out  to  a  radius  r=0.8325  5^  cm,  where  in  this  case  is  the 
initial  gradient  width.  This  value  of  r  is  the  Vi  downpoint  in  density. 
Therefore,  mass  and  energy  are  not  conserved  in  the  jet  system  due  to  shock 
heating.  Once  the  time  integration  is  complete,  the  simulation  then 
calculates  the  ratio  8rr/8rn  which  then  can  be  compared  with  results 
provided  by  pl2.c. 
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V.  SIMULATION  RESULTS 

A  INTRODUCTION 

Simulations  were  run  using  the  profiles  given  by  equations  (30)  and 
(31),  where 

8zr  =  1/3.95  cm  , 

8zn  =  1/4.0  cm  , 

and, 

n max(t )  —  1.2315xl0-3  r.u.  nsec-2 12  e~t/150  nsec  .  (34) 

The  maximum  density  function  (nmax(t))  was  obtained  by  fitting 
experimental  data  obtained  by  Bird.  The  experimental  data  of  nmax(t)  is 
shown  in  Figure  16.  The  modeled  mathematical  function  n audit)  is  shown  in 
Figure  19.  The  experimental  data  obtained  by  Bird  was  measured  at  r=0.4 
cm.  Since  this  is  approximately  the  e-1  downpoint  in  density  the  measured 
values  were  multiplied  by  three,  in  order  to  extrapolate  the  density  at  r=  0 
cm.  The  experimental  data  for  times  earlier  than  about  200  nsec  does  not 
accurately  represent  true  plasma  jet  conditions,  since  the  jet  has  not  arrived 
at  the  measuring  probe.  This  inaccuracy  was  neglected  in  the  initial 
simulations.  In  the  end,  modification  of  n  max(/ )  in  the  simulation  proved  to 
yield  no  new  information,  only  faster  magnetic  field  growth  rates.  All  the 
"useful"  physics  was  obtainable  using  the  mathematical  forms  given  above. 
The  characteristic  radial  widths,  8™  and  8ry,  were  varied  from  0.55  cm  to 
0.95  cm  in  order  to  observe  the  change  in  the  simulated  behavior.  Tmax  was 
also  investigated  to  find  an  estimate  of  the  post-shock  electron  temperatures 
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required  to  obtain  the  observed  magnetic  fields.  The  characteristic  widths  of 
the  axial  gradients  are  fixed  by  the  shock  heating  mechanism  to  be 
approximately  2-3  mm.  The  axial  temperature  gradient  width  (&zT)  is 
slightly  larger  than  the  axial  density  gradient  width  (82„)  due  to  the  high 
electron  heat  conduction  rates.  In  essence,  this  causes  electron  pre-heating 
in  front  of  the  shock.  Due  to  this  difference,  §2n  will  be  referred  to  as  the 
shock  thickness. 

Before  going  on,  it  should  be  noted  that  the  magnitude  of  the  shock 
generated  field  is  astonishingly  high  in  the  simulations.  The  simulations,  of 
course,  assume  no  losses,  when  in  fact  many  loss  mechanisms  exist.  This 
chapter  will  present  the  simulated  data  and  then  attempt  to  provide  some 
idea  of  the  loss  mechanisms  involved.  Comparison  with  experimental  data 
will  show  that  energy  losses  from  the  shock  heated  plasma  are  indeed 
significant. 

First,  the  effect  of  differences  in  radial  temperature  and  density 
diffusion  will  be  discussed.  It  will  be  shown  that  as  5rr/8m  increases,  the 
generated  magnetic  field  increases  and  the  axial  size  of  the  generated  field 
approach  those  observed  in  experiment.  Second,  the  background 
temperature  is  assumed  to  be  1.0  eV.  Since  the  plasma  jet  velocity  is  known 
from  experiment,  the  maximum  electron  temperature  due  to  shock  heating 
(Tmax)  is  43.6  eV.  This  high  temperature,  however,  is  above  the  dissociation 
and  ionization  potentials  of  N2,  N2,  N,  and  Al.  This  implies  that  only  a 
small  fraction  of  the  energy  produced  in  shock  heating  will  be  directed  to 
heating  electrons  and  creating  gradients.  Therefore,  section  C  will  examine 
how  electron  heating  efficiency  affects  the  shock  generated  magnetic  fields. 
Third,  the  temporal  behavior  is  examined.  It  will  be  shown  that  due  to  low 
electron  heating  efficiency  and  the  neglect  of  field  decay  due  to  collisional 
losses,  comparison  of  simulation  and  experiment  in  this  respect  is  not  good. 
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Finally,  the  field  geometries  are  compared  with  experiment  and  the  effects 
of  axial  electron  pre-heating  are  analyzed. 


Figure  19:  Mathematical  Model  for  nmax  vs.  Time  at  r=0  cm,  0  =  0°. 


B.  5rT  /  8rn  BEHAVIOR 

Figure  20  shows  the  results  of  studying  the  relationship  between  the 
maximum  generated  magnetic  field  at  400  nsec,  and  the  characteristic 
radial  temperature  and  density  widths.  As  usual,  positive  magnetic  fields 
represent  shock  generated  fields. 

The  abscissa  in  Figure  20  represents  the  ratio  of  radial  temperature 
width  to  radial  density  width.  There  are  several  features  to  notice  about 
Figure  20.  First,  if  8rT  =6m,  then  no  shock  generated  fields  are  produced. 
This  observation  immediately  supports  the  thesis  that  rapid  electron  heat 
conduction  in  the  radial  direction,  in  conjunction  with  shock  heating  at  the 
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Figure  20:  Generated  Magnetic  Field  versus  8^/8^.  £=400  nsec, 
8ot=0.25  cm,  S2r=l/3.95  cm,  Tbg=l.O  eV,  Tmax= 43.597  eV,  nmax=15  r.u., 
and  nbg  =  0.1  r.u..  (82y/82n  =1.01266). 


front,  produce  field  reversal.  Therefore,  drT  >  8m  is  a  fundamental 
requirement  in  achieving  field  reversal.  Clearly,  magnetic  field  generation  is 
driven  by  radial  heat  conduction  and  axial  shock  heating.  The  curve  in 
Figure  20  depicts  a  logarithmic  relationship  between  B  and  8^/8^,  which 
can  be  written, 


B 


8rT 


=  29345.7  Gauss  log  10 


8rT 


+  21.95  Gauss 


It  should  also  be  noted  that  the  magnetic  field  is  not  an  absolute  function  of 
the  ratio  8rT/8rn .  As  the  radial  size  of  the  plasma  jet  (§m)  increases,  for  a 
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given  8rr/5m ,  the  generated  magnetic  field  decreases  linearly. 

As  already  alluded  to  in  the  theory  chapter,  if  the  shock  forms  in  a 
background  with  no  initial  field  (laser  generated  field)  then  the  electrons  are 
not  magnetized.  The  collisionless  electrons  are  free  to  flow  away  from  the 
plasma  jet  to  form  an  electric  bi-layer  in  the  background  plasma.  This 
fieldless  mechanism  produces  an  initially  higher  8^/8™  and  hence  higher 
magnetic  field. 

Examination  of  the  experimental  data  shows  that  the  axial  width  of 
the  generated  magnetic  field  is  about  0.5  cm.  Figure  21  shows  that  when 
8 rr/8m  >1.0,  the  generated  field  indeed  approaches  a  constant  value  of  about 
0.46  ±  0.05  cm.  The  uncertainty  in  the  measurement  is  mainly  due  to  the 
subjectivity  in  measuring  the  field  width  from  the  computer  output.  When 
5rr/8m=1.0,  the  field  size  is  about  8.2  cm  then  as  radial  heat  conduction 
occurs  the  field  size  quickly  falls  to  observed  values.  This  observation 
supports  the  existence  of  differences  in  the  diffusion  rates  of  temperature 
and  density. 

Varying  maximum  electron  temperature  did  not  affect  the  reversal 
thickness.  8/t  and  together  with  the  degree  of  electron  pre-heating, 
predict  the  axial  dimensions  of  the  reversal!  This  is  remarkable  since  it  is 
not  immediately  obvious  that  the  axial  dimension  of  the  shock  generated 
reversed  field  is  independent  of  the  maximum  electron  temperature.  As  the 
ratio  brx  /  Srn  becomes  larger,  the  reversed  field  axial  width  essentially 
remains  constant  at  the  experimentally  observed  value. 

Simulation  pl7.c  was  used  to  determine  the  magnitude  of  the  ratio 
Srr/Sm.  The  following  model  was  used  in  pl7.c  to  determine  this  ratio.  As 
discussed  earlier,  the  electron  temperature  diffusion  is  modeled  by  Bohm 
diffusion,  which  is  a  model  that  describes  the  energy  transport  across  a 
magnetic  field.  As  calculated  earlier,  the  magnetic  field  is  the  dominant 
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mechanism  present  in  determining  electron  motion.  That  is,  the  electron 
collision  frequency  and  mean  free  path  are  are  small  enough  to  provide  gyro 
motion  of  the  electrons  about  the  magnetic  field  lines.  The  electrons  are 
inhibited  from  diffusing  freely  due  to  the  presence  of  a  magnetic  field  which 
modifies  their  motion.  It  is  important  to  understand  that  ions  are  not 
restricted  by  the  same  mechanism.  The  ions  have  a  much  lower 
gyrofrequency,  so  that  (ocxli  «  1;  hence,  they  are  less  inhibited  by  the 
presence  of  a  magnetic  field.  Therefore,  electrons  (due  to  their  lighter  mass) 
diffuse  their  thermal  energy  through  elastic  collisions  (in  gyromagnetic 
orbits)  with  other  electrons  at  a  very  high  rate.  Electron  heat  diffusion  is 
an  electron-electron  collision  dominated  mechanism  which  is  heavily 
dependent  on  the  presence  of  a  magnetic  field. 

Simulations  show  that  within  50  nanoseconds  8rr/Sm=1.13,  and 
continues  to  grow  for  the  remaining  duration  of  the  experiment.  Figures  9a 
and  9b  show  the  results  of  this  simulation.  One  essential  feature  of  pl7.c 
which  must  be  noted  is  that  between  r=0  and  the  Vi  downpoints  in  density 
and  temperature,  the  profiles  are  maintained  constant  through  the 
simulation.  This  is  indicated  in  Figures  9a  and  9b  by  vertical  lines  at 
approximately  r=0.4  cm.  This  was  done  because  in  reality  shock  heating  and 
compression  maintain  the  temperature  (Te (r ))  and  density  ( ne(r ))  within  the 
plasma  jet  as  long  as  the  jet  Mach  number  remains  constant,  which  is 
assumed  to  be  the  case.  Only  the  normal  component  of  velocity  is  shock 
heated,  hence  maximum  heating  and  compression  occurs  along  the 
centerline  r=0  cm.  The  temperature  and  density  within  the  plasma  jet  are 
functions  of  radius  due  to  the  curvature  of  the  shock  front.  If  the  shock  front 
holds  the  temperature  and  density  constant  to  a  radius  beyond  the  V2 
downpoint,  then  the  ratio  5ry/5r„  will  be  correspondingly  smaller. 
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The  numbers  ey  and  e„  in  Figures  9a  and  9b  give  an  indication  of  how 
much  energy  and  mass  must  be  added  to  the  plasma  jet  via  the  shock 
process  in  order  to  maintain  constant  profiles  of  temperature  and  density 
within  the  jet.  That  is,  by  50  nsec  energy  is  increased  by  15.1%  to  maintain 
constant  temperature  and  mass  is  increased  by  13.0%  to  maintain  constant 
density.  These  increases  must  be  made  to  offset  losses  by  diffusion. 

Simulations  show  that  the  ratio  5^/5^  increases  when  the  following 
parameters  increase;  maximum  electron  temperature  and  maximum 
electron  density.  The  ratio  8^/8™  decreases  when  the  background  electron 
temperature  and  density  increase. 

Given  that  Te>max=  43.6  eV,  7^  =1.0  eV,  /ie>max=15.0  r.u.,  and  n^=  0.1 
r.u.;  Figure  22  shows  the  growth  of  8^/5^  as  a  function  of  time.  At  0.5  nsec, 
8^/5^  =1.0.  Notice  that  the  slope  is  fairly  steep  up  to  about  100  nsec,  and 
then  decreases.  This  indicates  that  electron  temperature  diffusion  rate 
within  the  first  100  nsec  is  fast  and  then  slows  at  later  times.  This  may  be 
explained  by  the  fact  that  the  electrons  become  de-magnetized  outside  of 
r=l.l  cm.  Beyond  this  range  the  electrons  also  propagate  at  their  respective 
thermal  velocity.  At  700  nsec,  8^/8^  =1.31.  Again  these  values  depend  on 
the  range  to  which  the  shock  holds  the  temperature  and  density  constant 
within  the  jet. 

C.  Te  max  BEHAVIOR 

Maximum  electron  temperature  in  the  expanding  laser  plasma  is 
determined  by  the  Mach  number  of  the  shock.  Figure  2  shows  that  at  large 
Mach  numbers  (Af  x>2.4)  the  shock  heating  mechanism  becomes  significant. 
Maximum  electron  temperatures  on  the  order  of  40  eV  are  not  uncommon, 
depending  on  the  initial  background  electron  temperature.  pl2.c  was  used  to 
analyze  the  effect  of  maximum  electron  temperature  on  the  self-generated 
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Figure  21:  Axial  Thickness  of  Reversed  Shock-Generated  Magnetic 
Field  versus  8rT/5m.  For  8zT/8zn  =1.01266,  Szn=0.25  cm,  Tmax=43.6  eV, 
Tbg=  1  eV,  nmax=15  r.u.,  and  n^=0.1  r.u.. 

magnetic  field.  Figure  23  shows  the  result  of  investigating  the  response  of 
the  generated  magnetic  field  to  changes  in  electron  temperature.  The 
background  temperature  was  assumed  to  be  1.0  eV,  the  elapsed  integration 
time  was  400  nsec,  8rT  =  0.605  cm,  and  STO  =  0.55  cm  (8rr/S m  =  1.1).  8rr  was 
allowed  to  be  greater  than  §rn ,  so  that  the  affect  on  the  reversed  field  could 
be  studied.  As  will  be  discussed  in  the  next  section,  not  all  energy  in  the 
shock  heating  process  is  used  to  heat  electrons.  Ideally  (for  a  given  Mach 
number)  the  temperature  behind  the  shock  would  be  43.597  eV.  However, 
the  laser  and  background  plasmas  have  several  excitation,  ionization  and 
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Figure  22:  8^/5^  versus  Time  in  nanoseconds.  Assuming  Te  max  =43.6 
eV,  T^=1.0  eV,  ne>max=  15  r.u.,  and  n^=0.1  r.u.  This  simulated  data 
was  produced  by  pl7.c. 

dissociation  potentials  which  lie  below  this  maximum  theoretical 
temperature.  Therefore,  in  reality,  only  a  small  fraction  of  the  energy 
released  during  shock  heating  shows  up  as  an  electron  temperature 
increase.  This  inefficiency  reduces  the  observed  magnetic  field.  Figure  23 
shows  the  results  of  simulating  various  degrees  of  inefficient  electron 
heating.  The  abscissa  indicates  the  fraction  of  the  shock  heating  energy 
which  actually  gets  used  to  increase  the  electron  temperature.  The  ordinate 
indicates  the  fraction  of  the  magnetic  field  produced  relative  to  no  energy 
losses.  If  it  is  assumed  that  about  10%  of  the  energy  made  available  by 
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shock  heating  is  used  to  heat  electrons,  then  Figure  23  shows  that  the 
magnetic  field  generated  for  this  case  would  be  about  8%  of  the  total  no  loss 
magnetic  field.  Assuming  the  background  temperature  to  be  1.0  eV, 
maximum  temperature  to  be  43.597  eV,  background  density  to  be  0.2  r.u., 
maximum  density  to  be  15  r.u.,  5rj’/8rn=1.07,  and  =1.01266; 

simulations  show  that  the  magnetic  field  produced  was  910  Gauss. 
Accounting  for  92%  energy  losses  yields  82  Gauss.  It  should  be  noted  that 
this  analysis  does  not  account  for  the  ohmic  losses  in  the  plasma  which  are 
described  by  the  first  term  in  equation  (10).  The  maximum  observed  shock 
generated  field  in  experiment  at  this  background  pressure  (700  mtorr)  was 
about  60  Gauss  (see  Figures  15).  Therefore  (with  the  exception  of  ohmic 
losses),  90%  ionization,  excitation,  and  dissociation  losses  do  well  in 
reconciling  the  lossless  shock  magnetic  fields  simulated  by  pl2.c  with  those 
seen  in  experiment. 

It  should  also  be  noted  that  shock  field  generation  occurred  even  when 
the  generated  electron  temperature  corresponded  to  a  jet  Mach  number 
below  2.4.  This  result  leads  to  the  conclusion  that  strong  shock  fronts 
(Mi> 2.4)  are  not  necessary  in  creating  shock  field  phenomena.  However, 
large  Mach  numbers  do  provide  for  larger  shock  field  growth  rates.  In 
actuality,  resistivity  losses  cause  the  shock-generated  magnetic  fields  to 
decay  exponentially,  as  observed  by  McKee  [Ref.l]. 

D.  TEMPORAL  BEHAVIOR 

The  temporal  dependence  of  the  self-generated  magnetic  field  was 
studied  yielding  Figure  24.  Initially,  during  the  time  frame  when  the  laser  is 
still  irradiating  the  target  (t  ~  20  nsec),  SrT  ~  5 and  hzT  >•  82n .  Shortly  after 
laser  heating  has  begun,  however,  electron  heat  conduction  causes  the 
temperature  gradients  to  relax  faster  than  the  density  gradients.  Also 
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1 


Fraction  of  Energy  being  used  to  Heat  Electrons 


Figure  23:  Fraction  of  Total  Magnetic  Field  Produced  versus  Fraction 
of  Energy  being  used  for  Electron  Shock  Heating  (as  opposed  to  ioniza¬ 
tions  and  dissociations).  For  8^/8^  =1.1,  8zr/8zn  =1.01266,  8„,  =0.5477 
cm,  8zn=0.25  cm,  Tmax=43.6  eV,  7^  =1.0  eV,  nmax=15  r.u.,  nbg=  0.1  r.u.. 

within  50  nsec  after  cessation  of  the  laser  pulse,  a  shock  front  forms.  When 
this  occurs,  8rT>5rn  and  82r>82„ .  The  former  inequality  is  due  to  radial  heat 
conduction,  while  the  latter  inequality  is  due  to  axial  electron  pre-heating. 
As  this  condition  occurs,  shock  fields  are  generated.  Figure  24  shows  that 
the  magnetic  field  increases  linearly  with  time.  The  following  parameters 
were  used  to  derive  Figure  24:  8rr=0.55  cm,  8^/8^  =1.1,  8zn  =0.25  cm,  and 
8zT/&zn  =1.01266.  A  linear  fit  to  simulated  data  shows: 
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B(t)  =  3.6  Gauss  nsec  1 1  -  123.8  Gauss, 

where  t  is  measured  in  nanoseconds,  and  B  is  in  Gauss.  This  fit  is  valid 
between  50  nsec  and  700  nsec.  Of  course,  the  rate  of  magnetic  field 
production  is  dependent  on  5rr/8m ,  which  has  already  been  discussed;  it  is 
also  dependent  on  the  amount  of  electron  pre-heating  and  the  shock 
thickness.  This  dependence  will  be  discussed  in  section  F  of  this  chapter. 
Comparing  this  temporal  behavior  with  experiment  (Figure  24)  does  not 
show  a  good  agreement.  This  can  be  attributed  to  the  following  reasons. 
First,  the  simulation  only  accounts  for  field  production  by  temperature  and 
density  gradients  and  neglects  any  effects  due  to  field  diffusion  (ohmic 
losses).  Second,  not  all  of  the  energy  converted  from  kinetic  to  thermal  in 
the  shock  process  is  converted  to  an  increase  in  plasma  electron 
temperature.  Since  the  shock  temperature  is  about  43  eV,  excitation, 
ionization,  and  dissociation  of  the  nitrogen  background  and  the  plasma  jet 
occurs.  The  energy  used  for  ionizations  and  dissociations  is  not  available  for 
increases  in  electron  temperature;  thus  the  effective  increase  in  electron 
temperature  due  to  shock  heating  is  smaller  than  that  predicted  solely  by 
lossless  shock  theory.  Smaller  gradients  result  and,  hence,  smaller  fields  are 
generated.  Figure  25  shows  the  experimental  temporal  behavior  compiled 
from  Bird’s  data. [Ref. 2] 

Comparing  the  700  mtorr  background  pressure  to  5  mtorr  pressure 
shows  that  the  maximum  initial  magnetic  field  for  the  700  mtorr  case  (solid 
line)  decays  away  much  more  quickly  than  does  the  5  mtorr  case  (short 
dashed  line).  This  is  because  at  700  mtorr,  shock  heating  causes  reversal  to 
occur,  which  acts  to  "cancel"  the  initial  field  rapidly  when  compared  to  other 
decay  mechanisms  such  as  collisional  losses.  At  5  mtorr,  only  the 
conventional  loss  mechanisms  act,  so  that  field  decay  occurs  more  slowly. 
Figure  25  also  shows  that  shock  heating  begins  to  occur  at  a  time  between 
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Figure  24:  Shock-Generated  Magnetic  Field  versus  Time  (nsec).  (Simu¬ 
lation).  For  8zT/8zn  =1.01266,  8zn=0.25  cm,  8^/5^  =1.1,  8m  =0.5477  cm, 
T max=43.6  eV,  T^=1.0  eV,  nmax=15  r.u.,  n^=0.1  r.u.. 

30  nsec  and  100  nsec,  since  the  relative  peaks  at  about  50  nsec  are 
unaffected  by  this  mechanism.  Shock  generated  magnetic  fields  (in  700 
mtorr)  are  not  detected  until  300  nsec  (broad  dashed  line  B< 0  Gauss). 
Therefore  it  takes  approximately  200  nsec  to  algebraically  "cancel"  the  laser 
self-generated  magnetic  field.  Due  to  the  background  interactions  discussed 
earlier,  the  peak  magnetic  field  in  700  mtorr  is  expected  to  be  higher  than 
the  peak  field  in  5  mtorr. 
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Figure  25:  Self-Generated  Magnetic  Field  versus  Time  (nsec).  (Experi¬ 
ment)  Solid  Line  indicates  700  mtorr,  Short  Dashed  Line  indicates  5 
mtorr,  and  the  Broad  Dashed  Line  indicates  the  shock  generated  field 
at  700  mtorr.  Data  for  5  mtorr  is  extrapolated  beyond  300  nsec.  [Ref.2] 

E.  FIELD  GEOMETRIES 

The  evolution  of  the  shock  generated  magnetic  field  (simulated  by 
pl2.c)  is  shown  in  Figures  26a  through  26f.  These  runs  were  completed  at  50 
nsec  increments.  The  profiles  used  are  given  by  equations  (30)  and  (31), 
which  represent  a  background  pressure  of  700  mtorr.  The  characteristic 
gradient  widths,  drT  and  8m,  were  assumed  to  be  constant  at  5rT/5m=1.07, 
8,7,  =0.55  cm.  In  reality,  shock  heating  and  field  reversal  would  not  occur 
until  a  shock  has  formed.  This  occurs  at  some  time  between  30  nsec 
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(cessation  of  laser)  and  100  nsec  (first  detectable  evidence  of  field  reversal  in 
experiment).  The  times  represented  in  Figures  26  are  relative  to  when  the 
shock  formed,  not  when  the  laser  irradiation  began.  The  laser  generated 
magnetic  fields  (which  by  themselves  are  not  of  interest  in  this  thesis)  have 
been  neglected  in  the  simulation  (pl2.c)  since  diffusion  of  these  fields  for 
times  greater  than  30  nsec  could  not  be  taken  into  account.  Therefore, 
Figures  26  only  show  the  evolution  of  the  shock  generated  fields.  Once  a 
shock  has  formed,  8rT  >  8m ,  and  reversal  begins.  In  Figures  26a  through 
26f,  max  Te  denotes  maximum  electron  temperature  in  the  jet,  Bmax  and 
2?  min  indicate  the  maximum  and  minimum  magnetic  field  strengths  along  a 
line  at  r,  dt  is  the  numerical  time  increment  to  ensure  sufficient  time  steps 
are  taken  during  the  integration  process,  and  Z0  is  the  location  of  the  jet 
front.  The  upper  window  represents  a  two-dimensional  intensity  plot  of  the 
magnetic  field,  while  the  lower  window  represents  a  graph  of  magnetic  field 
strength  versus  axial  distance  along  a  line  at  r ,  where  r  is  given  at  the  right 
in  each  figure.  As  usual,  positive  field  values  indicate  shock  fields.  It  should 
also  be  mentioned  that  the  scales  of  the  abscissa  and  ordinate  for  each  of 
the  curves  in  Figures  26  change.  Therefore,  the  apparent  steepening  of  the 
magnetic  field  profiles  with  time  is  only  a  manifestation  of  the  scale  used. 
The  magnetic  field,  however,  does  exhibit  the  maximum  growth  rate  at  the 
plasma  jet  tip. 

Figures  27a  through  27e  show  experimentally  obtained  magnetic  field 
contours  at  250  mtorr  N2.  At  250  mtorr,  field  reversal  begins  at  about  300 
nsec.  At  700  nsec  the  reversed  field  has  completely  cancelled  the  initial  field 
component.  It  would  have  been  preferable  to  compare  experimental  data 
taken  at  700  mtorr  N2,  but  250  mtorr  was  the  highest  background  pressure 
where  B-field  contour  data  was  available. 
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Figure  26c:  Simulated  Magnetic  Field  Contours  at  200  nsec. 
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Figure  26d:  Simulated  Magnetic  Field  Contours  at  400  nsec. 
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Figure  26e:  Simulated  Magnetic  Field  Contours  at  500  nsec. 
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Figure  26f:  Simulated  Magnetic  Field  Contours  at  700  nsec. 
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Figure  27a:  Experimental  Magnetic  Field  Contours  at  0  nsec  in  250 
mtorr  N2  background.[Ref.l] 


Figure  27b:  Experimental  Magnetic  Field  Contours  at  40  nsec  in  250 
mtorr  N  2  background.  [Ref.  1] 
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Figure  27c:  Experimental  Magnetic  Field  Contours  at  120  nsec  in  250 
mtorr  N2  background.  [Ref.  1] 
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Figure  27d:  Experimental  Magnetic  Field  Contours  at  300  nsec  in  250 
mtorr  N2  background.  [Ref.  1] 
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Figure  27e:  Experimental  Magnetic  Field  Contours  at  700  nsec  in  250 
mtorr  N  2  background.  [Ref.  1] 


80 


Comparing  Figures  26  with  Figures  27  shows  general  agreement  in 
field  geometry,  however,  a  disparity  exists  in  field  strength  because  of 
differences  in  background  pressure  and  the  omission  of  decay  mechanisms 
in  the  simulations. 

F.  ELECTRON  PRE-HEATING  (5zr/5zn ) 

Electron  pre-heating  occurs  in  front  of  the  shock  due  to  electron  heat 
conduction  in  the  axial  direction.  As  the  magnitude  of  electron  pre-heating 
increases,  the  ratio  62r/8zn  of  course  also  increases.  This  change  affects  the 
magnitude  of  the  shock  fields  generated  and  the  size  of  the  field  itself. 
Figure  28  shows  that  as  the  degree  of  electron  pre-heating  increases,  the 
generated  magnetic  field  decreases.  This  figure  did  not  provide  any  insight 
into  the  actual  magnitude  of  5zr/52n  due  to  uncertainties  in  the  efficiency  of 
shock  heating.  However,  if  S^/S^  is  plotted  versus  the  field  size  (as  in 
Figure  29)  it  is  evident  that  because  experiment  showed  axial  field  width 
was  about  0.5  cm,  the  ratio  82y/5zn  must  be  less  than  about  1.05.  Therefore, 
5^/5^  =1.013  was  used  in  the  simulations.  Assuming  that  the  degree  of 
electron  pre-heating  is  small  (based  on  field  width  estimates)  and  that 
5ZT^zn  is  constant,  how  does  the  magnitude  of  the  generated  field  vary  with 
the  shock  front  thickness  (8ra)?  This  is  shown  in  Figure  30.  As  the  shock 
thickness  increases  by  50%,  the  generated  field  decreases  by  about  30%. 
Experiment  shows  that  the  shock  thickness  is  about  0.25  cm  to  0.3  cm. 
Therefore,  despite  the  unknown  magnitude  of  plasma  heating  losses,  this 
analysis  still  showed  that  the  magnitude  of  electron  pre-heating  is  small. 


81 


8  zT/8zn 


Figure  28:  Shock-Generated  Magnetic  Field  versus  the  Magnitude  of 
Electron  Pre-Heating  (S2r/82n).  For  8rr/8rn=l.l,  71max=43.6  eV,  Tbg=  1.0 
eV,  nmax=15  r.u.,  nbg=  0.1  r.u.. 
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Figure  29:  Axial  Reversed  Field  Thickness  versus  the  Magnitude  of 
Electron  Pre-Heating  (82r/S2„).  For  8^/8^  =1.1,  Tmax=43.6  eV,  7^  =1.0 
eV,  nmax=15  r.u.,  nbg=  0.1  r.u.. 
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Figure  30:  Shock-Generated  Field  versus  Shock  Front  Thickness  (5zn) 
for  a  constant  ratio  5zr/5zn.  For  5zr/8zfl  =1-01266,  5zn  =0.25  cm, 
8rr/8m=l.l,  5rn  =0.5477  cm,  Tmax=43.6  eV,  Thg-1.0  eV,  nmax= 15  r.u., 


VI.  APPLICATIONS 


A  INTERPLANETARY  SHOCKS 

The  purpose  of  this  chapter  is  to  provide  an  example  of  shock 
generated  magnetic  fields  which  may  be  observed  outside  the  laboratory. 
The  example  chosen  in  this  section  is  an  interplanetary  shock.  The  goal  will 
be  to  provide  an  estimate  of  the  steady  state  magnetic  fields  present  due  to 
shock  heating.  Because  the  specific  geometry  of  a  given  interplanetary  shock 
will  vary  greatly  and  collisionless  heat  conduction  mechanisms  are  involved, 
many  assumptions  will  have  to  be  made.  Detailed  analysis  of  the  shock 
generating  mechanism  in  a  collisionless  background  will  not  be  attempted 
here;  however,  this  would  provide  a  good  starting  point  for  further  research 
and  study.  The  analysis  performed  here  assumes  that  shock  generated  fields 
can  be  produced  in  a  collisionless  environment  since  the  existence  of 
collisionless  shocks  and  more  importantly  shock  heating  have  been  observed. 
Using  equation  (10)  and  neglecting  the  convective  term,  the  rate  of  change 
of  the  magnetic  field  will  be  assumed  to  be  zero.  This  will  be  the  case  when 
production  and  decay  mechanisms  achieve  an  equilibrium.  Therefore,  it 
becomes  possible  to  estimate  the  steady  state  magnetic  field  present  if  the 
gradients  and  conductivity  of  the  plasma  are  known.  First  order  estimates 
of  these  quantities  will  have  to  be  made. 

Large  fluctuations  in  the  activity  of  the  sun  induce  shock  conditions  in 
the  outward  propagating  solar  wind.  These  shock  waves  travel  outward  with 
the  solar  wind  at  velocities  approaching  800-1000  km/sec  when  they  leave 
the  sun’s  corona.  The  particles  present  in  the  solar  wind  are  mainly  protons 
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and  He**.  As  these  shocks  approach  the  earth,  satellites  measure  the 
discontinuity  in  the  plasma  density,  temperature,  and  velocity. 

Figure  31  shows  data  taken  by  Pioneer  11  in  May  of  1978  [Ref.13]. 
This  data  will  be  used  to  estimate  the  steady  state  value  of  the  expected 
shock  generated  magnetic  field.  Notice  that  the  temperature  shows  signs  of 
particle  "pre-heating"  prior  to  the  arrival  of  the  main  shock.  This  provides 
pivotal  clues  for  values  of  the  ratios  8rr/Sra  and  8zr/82n .  It  will  be  assumed 
that  this  "slight"  temperature  increase  can  be  attributed  to  axial  pre¬ 
heating,  as  discussed  in  the  previous  chapter.  The  first  objective  is  to 
determine  the  velocity  of  the  shock  so  that  times  can  be  related  to  distances. 
The  shocked  plasma  undergoes  a  velocity  increase  of  about  190  km/sec. 
Since  the  plasma  density  increased  by  a  factor  of  3.5,  this  implies  a  shock 
Mach  number  of  4.5  (assuming  y=1.66).  Now  the  velocity  of  ion  sound  in  the 
background  is  given  by: 

(y  +  1)M  x 
2 {Ml  -  1)  p 

where  up  is  the  "blast  wind"  velocity  behind  the  shock  which  was 
determined  to  be  190  km/sec.  Therefore,  c  x=59.4  km/sec  and  hence  the  shock 
velocity  is  266  km/sec.  When  referring  to  Figure  31,  one  hour  represents 
about  9.576x10s  km.  Various  gradient  widths  can  be  now  be  calculated: 
&zT  =  1.2xl06  km,  52n  =  5.7x10s  km.  This  implies  S^/S^  =2.09.  Since  no  data 
exists  for  the  radial  behavior,  it  will  also  be  assumed  that  =2.09. 

Letting  8„,  =  l.OxlO6  km,  it  is  now  possible  to  estimate  the  expected 
magnetic  field.  The  equation  to  be  solved  is: 


1  y2jg  +  k  dT  dn  _  dT  dtl  _  q 

|x0  a  n„e  dz  dr  dr  dz 


where  np  is  now  the  proton  density.  Substituting  the  partial  derivatives 
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with  approximate  expressions  involving  the  characteristic  gradient  widths 
gives: 


B 


H^oSrir  npe 


-(87V  8 nD) 


&rT  &zn  8m  hzT 


=  0. 


Let  87V  =  1.6x10®  K  and  8 nD  =  1  cm-3  =  10“®  m-3;  then  solve  for  B  yields: 


B(a)  =  (MooS^r) 


k 


npe 


(87V  5 nB) 


r 

1 

1 

8rr 

8m 

8m  8zr 


=  4.64xl0~2  GOm  o, 


where  B(o)  is  in  Gauss  and  a  is  in  mhos-m-1.  The  effective  conductivity  of 
the  plasma  (a)  is  still  unknown  and  is  difficult  to  calculate  for  the 
collisionless  interplanetary  example.  Figure  32  shows  a  plot  of  the  above 
relationship.  The  plasma  Spitzer  conductivity  runs  from  0  to  1000  mhos- 
m-1.  Therefore,  for  a  relatively  conductive  background,  a  steady  state 
magnetic  field  on  the  order  of  several  Gauss  is  expected.  If  the  background 
is  non-conductive,  then  the  magnetic  field  can  be  arbitrarily  small. 

In  conclusion,  although  the  specific  mechanisms  which  cause  the 
existing  gradients  in  a  collisionless  plasma  have  not  been  discussed,  it  is 
evident  that  the  general  phenomenon  which  occurs  in  the  laboratory  case  is 
also  applicable  in  the  interplanetary  case.  Therefore,  shock  generated 
magnetic  fields  in  astrophysical  phenomena  can  be  measured  if  the 
background  is  sufficiently  "conductive". 
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Figure  31:  Interplanetary  Shock  data  taken  from  Pioneer  11  on  11  May 
1978  [Ref.  12]. 
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Plasma  Conductivity  (mhos/m) 


Figure  32:  Shock  Generated  Magnetic  Field  versus  Background  Plasma 
Conductivity  for  an  Interplanetary  Shock. 

B.  NUCLEAR  EMP  EFFECTS  IN  THE  MHD  DOMAIN 

High  altitude  (>100,000  ft)  explosions  of  nuclear  weapons  produce 
electromagnetic  signals.  One  such  signal  is  called  the  electromagnetic  pulse 
(EMP).  A  typical  EMP  waveform  of  a  high  altitude  nuclear  burst  is  shown  in 
Figure  33.  Although  actual  EMP  waveforms  are  classified,  Figure  33 
provides  accurate  order  of  magnitude  data.  There  are  basically  two  regions, 
the  prompt-gamma  signal  and  MHD  signal.  For  times  less  than  1  second, 
prompt  gamma  and  neutron  scattering  produce  a  large  well  defined  pulse. 
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After  1  second,  late-time  EMP  is  evident  which  is  produced  by  MHD 
phenomena.  Since  the  discovery  of  EMP  in  the  early  1950’s  during  nuclear 
tests,  much  theoretical  and  experimental  work  has  been  completed  on 
understanding  the  prompt-gamma  signal.  This  early  signal  can  cause 
significant  damage  to  electronic  equipment  and  therefore  is  of  extreme 
military  interest.  The  late-time  MHD  (Magnetohydrodynamic)  signal, 
however,  has  not  been  so  widely  studied  and  is  still  not  fully  understood. 
Although  immediate  equipment  damage  by  this  signal  may  not  occur, 
certain  systems  such  as  long  cable  or  wiring  systems  (submarine  ELF 
antenna,  telephone,  etc...)  may  shutdown  due  to  overload  protection  devices. 
Large  potentials  on  the  order  of  several  kilovolts  per  100  km  may  result. 
The  focus  of  this  section  will  be  to  show  that  the  electromagnetic  signature 
associated  with  a  nuclear  EMP  in  the  MHD  domain  can  be  explained  in 
terms  of  shock  generated  fields!  The  success  of  the  following  analysis  is  a 
testament  to  the  fundamental  nature  of  the  shock  field  mechanism  and  to 
its  great  wide  ranging  application.  The  late-time  generation  of  an 
appreciable  electric  field  is  the  direct  focus  of  this  section. 

It  is  assumed  that  a  250  kiloton  nuclear  warhead  is  detonated  at  an 
altitude  of  100,000  ft  (30,500  m).  At  this  altitude  only  about  60%  to  90%, 
depending  on  atmospheric  conditions,  of  the  yield  is  transferred  to  shock 
energy.  The  remaining  energy  is  coupled  into  the  initial  release  of  thermal 
radiation  and  soft  x-rays.  As  the  altitude  increases,  coupling  of  the  energy 
yield  into  shock  energy  decreases  because  the  atmospheric  density  decreases 
with  increasing  altitude.  Within  tens  of  microseconds  after  detonation,  a 
intense  flux  of  x-rays  and  gamma  rays  fully  ionizes  the  surrounding 
atmosphere  out  to  several  kilometers.  Based  on  atmospheric  densities  at 
this  altitude,  the  electron  particle  density  is  4.7xl023  m-3  [Ref.  14].  Compton 
and  photocurrents  generate  the  short-time  EMP  effects  mentioned  earlier. 
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The  electrons  and  surrounding  ions  have  sufficient  time  to  achieve  thermal 
equilibrium  before  the  dynamic  shock  wave  arrives.  It  will  be  assumed  that 
the  surrounding  atmosphere  is  isotropically  heated  by  the  initial  thermal 
radiation  and  that  the  surrounding  temperature  is  increased  by  Tc. 
Therefore,  after  several  milliseconds  the  atmosphere  assumes  a  temperature 
of  Tc  +Tamb,  where  Tamb  is  the  local  ambient  temperature.  Tamb  (r,0)  is  a 
function  of  position.  That  is,  it  depends  on  the  height  above  the  earth’s 
surface. 

In  order  to  perform  calculations,  spherical  coordinates  will  be  used. 
The  coordinate  system  is  centered  about  the  burst  location,  r  indicates  the 
range  from  the  burst  point,  0  measures  the  polar  angle  from  the  zenith,  and 
<(>  measures  the  azimuthal  angle  around  the  horizontal  plane.  It  is 
immediately  apparent  that  the  problem  indicates  azimuthal  symmetry  and 
that  the  electron  temperature  and  density  will  not  be  dependent  on  <j>. 
Figure  34  shows  the  geometry  used  and  the  gradients  which  exist  in  this 
problem. 

According  to  reference  14,  the  shock  wave  of  a  250  KT  burst  achieves 
a  1  to  2  mile  radius  in  approximately  5  seconds.  It  will  be  assumed  that  the 
expanding  spherical  shock  obtains  a  radius  of  1.5  miles  (2414  m)  in  5 
seconds.  Therefore,  the  average  shock  velocity,  vsh ,  is  1584  ft/s  (482.8  m/s). 
At  5  seconds  the  shock  has  a  diameter  of  3  miles  (4828  m).  Due  to  the 
exponential  thinning  of  the  earth’s  atmosphere  as  a  function  of  height,  the 
temperature  and  density  at  a  height  of  2414  m  above  the  burst  point  are 
different  than  at  a  height  of  2414  m  below  the  burst  point.  Accepted  values 
of  ambient  temperature  and  density  2414  m  above  the  burst  point  are;  231 
K  and  0.0153  Kg/m 3,  respectively.  At  2414  m  below  the  burst  point  the 
ambient  temperature  and  density  are;  226  K  and  0.0244  Kg/m3,  respectively. 
It  should  be  noted  that  the  vertical  density  and  temperature  gradients  are 
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oppositely  oriented  and  that  the  density  gradient  is  greater  than  the 
temperature  gradient.  After  radiation  heating  has  occurred,  the  upper  and 
lower  temperatures  become,  (Tc  +  231  K)  and  (Tc  +  226  K),  respectively. 


Since  sound  speed  is  a  function  of  temperature,  the  shock  Mach  number  will 
also  be  a  function  of  temperature.  The  Mach  number  above  the  burst  point 


(2414  m)  will  be: 

Mt  = 


Vsh 

231  K 

cs,T 

Tc  +  231  K 
.  - 

Yt 


=  1.583 


231  K 


Tc  +  231  K 

and  the  Mach  number  below  the  burst  point  will  be: 


Vi 


Mb  = 


Vsh 

226  K 

CSJi 

Tc  +  226  K 

Vi 


=  1.601 


226  K 


Tc  +  226  K 


Vi 


where  cStT  (=304.9  m/s)  and  cSiB  (=301.6  m/s)  are  the  ambient  sound  speeds 
for  the  temperatures  231  K  and  226  K  above  and  below  the  burst  point, 
respectively.  The  Mach  number  below  the  burst  point  is  larger  than  the 
Mach  number  above  the  burst  point  simply  because  the  ambient  sound 
speeds  are  different.  The  shock  velocity,  however,  is  assumed  to  be  isotropic. 
The  vertical  temperature  gradient  in  the  atmosphere  implies  asymmetric 
shock  heating.  It  should  be  noted  that  in  order  to  maintain  M>  1,  the 
maximum  increase  in  temperature  produced  by  radiative  heating  is  about 
347  K.  This  implies  that  the  unshocked  background  can  have  a  maximum 
temperature  of  578  K.  Using  equation  20  to  now  calculate  the  shock  heated 
temperatures  (using  y  =  cp  /cv  =  1.4)  gives: 

T2>t  =  -2.4x1 0-4  K"1  T?  +  0.875TC  +  318.1  K 


and 


T2J}  =  -2.4X10-4  K"1  Tc2  +  0.836TC  +  313.3  K 
where  T2j  and  T2$  are  the  shock  heated  temperatures  above  and  below  the 
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burst  point,  respectively.  See  Figure  34.  Furthermore,  the  temperature 
difference  between  the  two  pole  locations  is: 

AT  —  T 2P  ~  T 2j}  “  4.8  K  —  O.OOIT*  . 

The  difference  between  the  shocked  and  unshocked  temperature  ( 8Te )  is 
then: 

8Te  =  T%t  -T1>t  =  -2.4X10"4  K"1  Tc2  -  0.165TC  +  87.1  K  , 

where  T1t=Tc  +  231  K.  It  will  be  assumed  that  6 Te  is  isotropic,  although  in 
reality  it  may  vary  slightly  as  a  function  of  altitude. 

Similarly,  using  the  relation  between  the  shocked  and  unshocked 
density  given  by  equation  18,  the  shocked  density  above  the  burst  becomes 
(assuming  the  unshocked  density  (nir)  is  3.6xl023  m-3  [Ref.14]): 

„  _  5xl026  K/m~3 

2>T  2 Tc  +  693.5  K  ’ 

where  u2it  is  in  m-3  and  Tc  is  in  Kelvin.  Furthermore,  letting  the 
unshocked  density  below  the  burst  (n  13)  be  5.7X1023  m~3  the  shocked 
density  becomes: 

=  7.98X1026  K/m~3 
23  2 Tc  +  683.6  K  ’ 


Finally,  the  vertical  density  gradient  due  to  atmospheric  thinning  is: 
An  =  n2tT  ~  n2$  -  -5xl026  K/m-3 


1.192TC  +  423.2  K 


(2 Tc  +  693.5  KX2TC  +  683.6  K) 
and  the  difference  in  shocked  and  unshocked  density  (6 ne )  is: 


3 ne  =  n-2tT  —  nip  = 


3.6xlQ23  m~3 
2 Tc  +  693.5  K 


(695.4  K  -  2TC ) . 


One  additional  term  remains  to  be  calculated  in  order  to  estimate  the 
magnetic  field  generated.  This  is  the  conductivity  of  the  atmospheric  plasma 
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(a).  Using  the  Spitzer  resistivity,  r|,  gives  [Ref.15]: 


r\  1.03X10-4  Q  m  eV3/2  Z  In  A 

where  T  is  in  eV,  Z  is  the  ionization  state,  and  In  A  is  the  coulomb 
logarithm.  Converting  eV  to  Kelvin  and  letting  Z-\  and  lnA=10  in  the  above 
equation  gives: 

o  =  7.77XKT4  mhos  m"1  eV"3'2  Tm  , 


where  o  is  in  mhos/m.  In  this  specific  case,  T  -Tc  +  228  K,  where  228  K  is 
the  ambient  temperature  at  the  burst  height. 


The  magnetic  field  is  now  calculated  by  a  method  similar  to  that  used 
in  the  previous  section.  The  convection  term  in  the  magnetic  field  generation 
equation  (eq.10)  is  assumed  to  be  negligible.  It  is  desired  to  solve  for  the 
magnetic  field  when  its  time  rate  of  change  is  zero,  i.e.  the  maximum  field. 
Therefore,  using  spherical  coordinates,  the  field  equation  becomes: 


dg_ 

dt 


=  0  =  — v2b  + 


* 

k  1  dT  dn  _  dT  dn 
ene  r  dr  30  30  3r 


It  is  difficult  to  work  in  terms  of  the  coordinate  0  at  this  point  in  the 
calculation.  Therefore,  a  change  to  the  altitude  variable,  h,  will  be  made. 
The  variable,  h,  represents  the  vertical  height  above  or  below  the  burst 
point.  The  relation  between  h  and  0  is: 


h  =  rcos0  , 


where  r=1.5  miles  (2414  m)  in  this  case.  Now  the  partial  derivative  with 
respect  to  0  can  be  written: 


_3_ 

30 


=  -r  sin0 


_3_ 

3 h  ’ 


Inserting  this  back  into  the  field  equation  gives: 


94 


o  =  _1_v2b  - 


enP 


-sin0 


dT  dn  dT  d n 


dr  dh  dh  dr 


Now  approximating  the  partial  derivatives  as  gradient  ratios  (as  done  in  the 
previous  section)  gives: 


0  = 


1  B 
M-oCT  5| 


k  sin0 

STe  a n 

AT 

ene 

8rj’  Ah 

* 

Ah 

where  8B  is  the  characteristic  size  of  the  magnetic  field  (8S  =  M/2),  M  is  the 
characteristic  width  of  the  vertical  density  and  temperature  gradient  of  the 
atmosphere  (4828  m),  8Te  is  the  shock-unshocked  temperature  difference 
calculated  earlier,  An  and  A T  are  the  vertical  density  and  temperature 
change  experienced  in  a  height  change  Ah  about  the  burst  point  after  it  has 
been  shock  heated,  and  finally  5ne  is  the  change  in  density  across  the  shock 
front.  Note  that  the  factor  of  sin0  ensures  that  magnetic  flux  is  only 
generated  where  the  temperature  and  density  gradients  are  perpendicular. 
At  0=0  and  0=n,  the  atmospheric  and  shock  gradients  are  parallel  and  hence 
do  not  contribute  to  the  field  production.  However,  around  the  "equator"  of 
the  expanding  shock,  the  gradients  are  perpendicular  and  result  in 
maximum  field  generation.  The  magnetic  fields  are  azimuthal  (-$  direction), 
due  to  the  symmetry  in  that  direction. 

Solving  for  the  magnetic  field  in  the  above  equation  and  assuming 
0  =  jt/2  gives: 


hjx0 a  Ah 
4 ene  8rT 


Sne  AT  -  8 Te  An 


(34) 


If  Tc  =340  K,  then  all  over  parameters  are:  ne  =  4.8xl023  m-3,  M=4828  m, 
5rT  =  1  m,  8ne  = -4xl021  m~3,  An  = -2.2xl023  m-3,  AT=4.5  K,  o=10.5  mhos/m, 
and  871e=-3.3  K.  This  all  yields: 


B  =  -2.08xl0"2 


1  -  0.025 


Gauss  . 
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B  is  in  Gauss.  Notice  that  the  behavior  of  this  equation  differs  significantly 
from  the  behavior  displayed  by  any  previous  analysis  of  shock  generated 
fields.  Shock  generated  fields  are  produced  even  when  5^/5^  =1.  This  was 
not  the  case  for  laser  produced  plasma  jets.  This  can  be  explained  by  the 
fact  that,  in  the  case  of  the  nuclear  burst,  the  vertical  atmospheric  gradients 
are  independent  of  the  radial  shock  gradients  generated  by  the  explosion. 
Non-parallel  gradients  exist  regardless  of  whether  the  electron  temperature 
propagates  outward  faster  than  the  electron  density.  Although  the  "bi- 
diffusive"  phenomenon  may  still  operate  in  this  regime,  it  is  not  vital  to 
generating  the  magnetic  field. 

Figure  34  shows  a  plot  of  EMP  generated  field  versus  initial 
background  heating  by  the  thermal  radiation  (Tc).  It  has  already  been 
shown  that  radiative  heating  can  not  exceed  about  340  K.  This  implies  that 
the  pre-shock  heated  atmospheric  temperature  can  not  be  more  than  about 
570  K.  Experimental  observation  shows  that  the  shock  wave  does  not 
dissipate  until  long  after  5  seconds  have  passed.  Figure  35  is  the  graph  of 
equation  34  when  all  terms  and  factors  involving  Tc  are  used.  The  full 
equation  is  lengthy.  Writing  it  down  would  not  add  anything  constructive  to 
this  analysis.  This  graph  shows  that  the  maximum  generated  field  will  be 
about  0.22  Gauss  when  Tc=150  K  (Tamb=380  K).  This  graph  also  assumes 
that  the  ratio  5rT/8m  =1. 

The  maximum  generated  field  is  about  0.22  Gauss.  Since  this  field  is 
time  varying  and  is  generated  in  a  period  of  about  5  seconds,  an  electric 
field  also  exists.  Using  Maxwell’s  equation,  the  electric  field  may  be 
calculated.  The  equation  of  interest  is: 

Vx£f  =  -  4^- . 
dt 

Approximating  this  equation  gives: 
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where  8#  is  the  characteristic  width  of  the  electric  field  (dE  =  8B  =  Ah  12). 
Using  AS  =0.22  Gauss,  A*  =5  seconds,  and  8#  =2414  m,  yields  E=  10-2  V/m. 
Furthermore,  because  the  magnetic  field  is  in  the  direction,  the  cross 
product  implies  the  electric  field  is  in  the  f  and  G  directions.  These  are  the 
proper  field  directions  observed  in  tests.  Comparison  of  this  result  to  actual 
field  values  measured  during  tests  shows  that  they  are  in  agreement. 
Although  this  analysis  has  been  rather  "crude"  in  that  many  assumptions 
were  made,  none  of  them  have  been  unreasonable.  The  exact  degree  of 
radiative  heating  is  not  known;  however,  simple  estimates  do  indicate  that 
appreciable  fields  can  be  generated  by  the  gradients  described. 

This  chapter  shows  that  interplanetary  shocks  and  nuclear  EMP 
effects  generate  magnetic  fields  using  a  similar  process.  This  conclusion  is 
quite  remarkable.  Further  study,  and  research  into  additional  applications, 
is  warranted.  A  few  additional  applications  are  given  in  the  concluding 
chapter  of  this  thesis. 
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Figure  33:  A  typical  High  Altitude  EMP  Waveform.  The  Electric  Field 
in  V/m  versus  Time  after  burst  in  seconds.  Notice  scale  is  log-log.  The 
maximum  Electric  field  in  the  MHD  domain  is  about  0.01  V/m. 
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Figure  34:  Geometry  of  a  High  Altitude  EMP  Burst.  Dominant 
gradients  are  depicted  in  the  figure.  The  shock  radius  is  1.5  miles 
(4824  m)  at  5  seconds.  [Ref.  14] 
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Figure  35:  Shock  Generated  EMP  Field  as  a  function  of  Atmospheric 
Radiative  Pre-heating.  {Big)  vs  Tc)  The  magnetic  field  is  an  extrema  at 
about  Tc  =150  K.  This  implies  unshocked  atmospheric  temperatures  are 
about  380  K. 
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VII.  CONCLUSION 


Shock-generated  magnetic  fields  produced  by  a  plasma  jet  propagating 
in  an  ionized  background  have  been  simulated  using  temperature  and 
density  profiles  derived  from  experimental  data.  Simulation  results  lead  to 
several  conclusions.  First,  shock  heating  occurs  at  the  jet  front  at 
appreciable  background  pressures  (>200  mtorr).  The  jet  has  a  velocity  of 
about  Mach  11.75  in  a  background  pressure  of  700  mtorr.  At  these  high 
Mach  numbers,  the  relative  temperature  rise  across  the  shock  front  exceeds 
the  density  increase  (see  Figure  2).  Simulations  show  that  the  shock 
mechanism  can  create  magnetic  fields  at  the  plasma  jet  front  even  at  low 
Mach  numbers  (Afx  <  2.4).  In  fact,  it  seems  that  there  are  two  criteria  for 
creation  of  shock  magnetic  fields:  (1)  Supersonic  plasma  flow,  i.e.  (Mi  >  1.0), 
and  (2)  a  difference  between  radial  temperature  and  density  diffusion.  Prior 
to  these  investigations,  electron  temperatures  were  largely  unknown. 
Assuming  a  background  of  1  eV  implies  maximum  electron  temperatures 
are  on  the  order  of  43  eV.  Simulations  confirm  that  magnetic  fields 
generated  by  shock  heating  are  azimuthal  with  respect  to  the  jet.  At  larger 
background  pressures  (>  200  mtorr),  field  reversal  occurs  at  the  jet  front 
when  radial  electron  heat  conduction  occurs  faster  than  radial  ion  density 
diffusion.  Density  diffusion  is  dictated  by  ion  motion  due  to  the  larger  ion 
mass  and  the  requirement  of  quasi-neutrality.  The  directed  motion  of  the 
ions  in  the  axial  direction  is  much  larger  than  their  radial  thermal  motion. 
Even  though  the  electrons  are  magnetized,  the  electron-electron  interactions 
provide  that  radial  thermal  transport  occurs  more  rapidly  than  ion 
transport.  This  bi-diffusivity  of  the  electron  density  and  temperature  radial 
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gradients  causes  field  generation  in  the  reverse  direction  at  the  tip  of  the 
jet.  In  addition  to  radial  heat  conduction,  axial  heat  conduction  also  occurs 
which  pre-heats  electrons  in  front  of  the  shock.  This  pre-heating  is  shown  to 
be  small  where  8zT/8zn  £  1.05.  Approximating  the  axial  and  radial  gradients 
by: 


v2re 

bzT 


V2ne 


ne 


V  T 

y  r  ±  e 


8rT 


Vr^e  = 


shows  how  the  reverse  magnetic  field  production  due  to  shock  heating  is 
accomplished.  The  requirement  to  produce  shock  generated  magnetic  fields 
is: 

IV2Te  I  IVrne  I  _  (5^X5^) 
lVrTel  IV2n,l  (8zT)(8m)  • 

The  5s  are  the  gradient  scale  lengths.  Shock  heating  at  the  jet  front  implies: 

8zt  ~  8zn  and  8zx  <  5ry  . 


Therefore  electron  heat  conduction  implies: 

5rr  >  5m  . 

Simulations  confirm  this  result.  Simulations  predict  that  axial 
dimension  of  the  shock  generated  magnetic  field  is  only  a  function  of  8zT  and 
5^.  That  is,  if  8rT  /  5rn>1.0,  then  as  8zT/8zn  increases  the  field  width  also 
increases.  This  relationship  is  shown  in  Figure  29.  Experimentally,  axial 
reversal  thickness  was  about  0.5  cm  which  corresponds  to  the  condition: 

5 rT  >  8rn  and  52y  <  1.05  52n  . 

Maximum  electron  temperature  and,  hence,  jet  Mach  number  do  not 
affect  the  shock  field  width  (reversal  thickness).  Experimentally,  shock 
heating  occurs  within  100  nsec  after  the  cessation  of  the  laser  pulse.  Once 
shock  heating  occurs,  the  growth  of  the  shock  generated  magnetic  field 
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occurs.  For  background  pressures  less  than  200  mtorr,  shock  formation  does 
not  occur  (or  occurs  too  late)  and  free  expansion  of  the  jet  prevails  (see 
Figure  12).  The  simulation  program  (pl2.c)  neglects  the  convection  and 
diffusion  terms  in  equation  (10),  as  well  as  the  ionizations,  excitations,  and 
dissociations  which  occur.  Therefore,  with  the  exception  of  the  temporal 
dependence  and  magnitude  of  the  magnetic  field,  pl2.c  provides  good 
agreement  with  experiment. 

In  this  instance,  numerical  simulations  (guided  by  experimental  data) 
provided  useful  verification  and  insight  into  the  shock  heating  mechanism  in 
plasmas.  pl2.c  also  provided  an  increased  basic  understanding  as  to  which 
jet  and  background  parameters  were  relevant  and  how  each  affected  the 
mechanism  of  shock  generated  magnetic  fields.  Applicability  of  this  basic 
study  is  wide  ranging;  some  areas  might  include:  (1)  explaining  late-time 
electromagnetic  pulse  (EMP)  effects  after  a  nuclear  detonation,  (2)  magnetic 
field  phenomena  associated  with  jets  created  by  Rayleigh-Taylor  instabilities 
(inertial  confinement  fusion),  (3)  streaming  jets  found  in  astrophysical 
phenomena,  (4)  magnetic  fields  associated  with  interplanetary  shocks 
formed  by  the  solar  winds,  (5)  the  earth’s  bow  shock,  (6)  the  dynamics  of 
comet  tail  formation  and,  (7)  electromagnetic  signature  of  an  orbital  re-entry 
vehicle.  Further  study  in  these  fields  in  relation  to  the  shock  mechanism 
seems  warranted. 

Two  specific  applications  were  analyzed  in  this  thesis;  interplanetary 
shocks  and  late-time  nuclear  EMP  fields.  The  shock  fields  generated  by  the 
interplanetary  shock  vary  depending  on  the  conductivity  of  the 
interplanetary  plasma.  Although  the  second  application  seems  completely 
unrelated  to  the  first,  the  same  field  production  mechanism  fundamentally 
links  both  phenomena.  Late-time  EMP  fields  are  also  generated  by  non¬ 
parallel  gradients.  Calculations  show  that  vertical  atmospheric  gradients 
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crossed  with  the  radial  shock  gradients  produce  magnetic  and  electric  fields 
which  agree  well  with  those  measured  during  tests.  The  fact  that  such 
widely  varying  phenomena  are  described  by  a  single  mechanism  is  a 
testament  to  the  predictive  power  and  fundamental  nature  of  this  physical 
process. 

Magnetic  field  generation  by  the  dynamo  mechanism  has  been  studied 
extensively  in  the  terrestrial  and  astrophysical  domain;  however,  a  second 
mechanism  is  also  prevalent  in  the  universe.  The  shock  generation 
mechanism  described  in  this  thesis  may  help  to  explain  many  examples  of 
dynamic  magnetic  field  production  in  the  universe.  In  particular,  those 
situations  where  gas  or  plasma  jets  stream  supersonically,  so  that  shock 
heating  produces  non-parallel  temperature  and  density  gradients.  An 
example  of  this  has  been  thoroughly  studied  in  the  laboratory.  Nature 
displays  extravagant  phenomena  of  this  sort  as  well. 
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APPENDIX 


A.  PL2.C 


/*★★★★*★★★★*******★★★*•***★*★★★★★★★★★★**★**★*★★*★★★ 
====  Self  Generated  Magnetic  Fields  ==== 
PL2.C 

tinclude  <stdio.h> 
tinclude  <cmath.h> 
finclude  <greeks.c> 

ARRAY 2d  B,  T,  n; 


float  t,z,r,z0,a, nmax, dz, dr,b,  dt, rl, tmax, c [ 40] , Bgrmx, Bgrmn, 

Bgrx,  dtdr, dtdz, dndr, dndz, tl,bmn,bmx,DZt,DZn, Tmx,Tbg, 
aT, an; 

int  zck, nc,  i,  nz,  nr,  j, Bint, nrl,  it; 
main()  { 

/*==================  read  in  parameters  ================—=*/ 


printf ("\n***  B-field  intensities  ***\n") ; 

printf("  r,  tmax,  nc,  DZt,  DZn,  Tmx,  Tbg,  aT,  an\n  ")  ; 

scanf ("%f %f%d%f%f %f %f %f %f ", &rl, &tmax, &nc, &DZt, &DZn, &Tmx, &Tbg, &aT, San) ; 

tl=0 . 0062*tmax; dz=t 1/250 . 0; dr=3 .0/250 . 0;dt=0 . 02; z 0=0 . 0031*tmax; 

DZn=l/DZn;DZt=l/DZt; 

/*====================  Initialize  Arrays  ==================*/ 

for  (nz=0;nz<250;++nz)  {  for  (nr=0;nr<250;++nr)  { 

B[nz] [nr]=T[nz] [nr]=n[nz] [nr]=0.0;}} 

for  (t=0 . 0; tctmax; t=t+dt)  { 

zck- (int) ((250. 0/tl)*((-0. 0031 *t)+z0) ) ; 

nmax=0 . 0012315*t*t*exp (-t/150 . 0) ; 

/*===========  Initialize  T  and  n  arrays  for  time  t  ========*/ 

for  (nz=0;nz<250;++nz)  {z=(tl/250 .0) * (float)nz; 

for  (nr=0;nr<250;++nr)  { 

r- (1. 5/125. 0).  *  (float  )nr-l.  5; 

if  (z>=z0){  T[nz] [nr] =Tbg;n [nz] [nr]=0.2; } 

if  (z<z0) 

{T [nz] [nr]= ( (Tmx-Tbg) * (1 . 0-exp (aT* (z-zO) ) ) *exp (-r*r*DZt*DZt ) ) +Tbg; 
n [nz]  [nr]  =  (nmax* (1 . 0-exp (an* (z-zO) ) ) *exp (-r*r*DZn*DZn) ) +0 .2; } } } 
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/*==============  find  derivatives  wrt  z  and  r  ==============*/ 

for  (nz=l;nz<249;++nz)  { 

if  (nz>=zck)  { 

for  (nr=l;nr<249;++nr)  { 

dtdz=  (T  [nz+1]  [nr]  -T  [iiz-1]  [nr] )  /  (2 . 0*dz) ; 

dndz= (n [nz+1] [nr] -n [nz-1] [nr] ) / (2 . 0*dz) ; 

dtdr= (T [nz] [nr+1] -T [nz] [nr-1] ) / (2 .0*dr) ; 

dndr= (n [nz] [nr+1] -n [nz] [nr-1] ) / (2 . 0*dr) ; 

/ *================  Calculate  updated  B-field  ===============*/ 

b=l/n [nz] [nr] ; 

B[nz] [nr]=B[nz] [hr] +(b* (dtdz*dndr-dtdr*dndz) ) *dt; } } } } 


/*=================  Find  Bgr  max/min  =======================*/ 

Bgrmx=Bgrmn=0 . 0 ; 

for  (i=l; i<249; ++i)  {  for  ( j=l; j<2 4 9 ; ++  j )  { 
if  (Bgrmx<=B[i]  [ j] )  Bgrmx=B [i] [ j] ; 
if  (Bgrmn>B [ i ]  [ j ] )  Bgrmn=B [i]  [ j ] ;  } } 
if  (Bgrmx>=-Bgrmn)  Bgrx=Bgrmx; 
if  (Bgrmx<-Bgrmn)  Bgrx=-Bgrmn; 

init();  color_scale  ( "bluered" ) ;  grey_scale  ("greyscalel")  ;  windowOO; 
bgcol (7) ;  erase  ();  color (0);  mode("Fl"); 
rect (237, 80, 841,195) ;  color(2); 

rect (230, 75, 834,190) ;  mode("F0");  color(O);  rect (230, 75, 834, 190) ; 
move (310, 80) ;  color (4); 

printf ("**  Magnetic  field  Intensity  **");  color (0); 
move  (280, 115) ;  printf  (,,Time=%g  nsec",tmax); 
move(280, 150) ;  printf ("P") ;sub('a' ) /printf ("=700  mtorr") ; 
move (110, 200) ;  printf  ("1.5") ;  move (85, 325) ;  printf ("r (cm)  ") ; 
move (95, 463) ;  printf ("-1. 5") ; 

move (170, 720) ;  printf  ("0") ;  move (395,  720)  ;  printf("z(cm)"); 
move(635,  720)  ;  printf  ("%2  .  If  ",tl)  /  windowOO; 

/*===========  b-z  graphics,  B-field  intensity  ===============*/ 

if  (nc ! =0 )  {for  (i=0; i<nc/++i)  c [i] =Bgrx* (i+0 . 5) / (f loat ) nc; 
for  (i=0; i<250; ++i)  for  ( j=0; j<250; ++ j )  if  (B [i]  [ j ] <0 . 0 ) 

B [i]  [ j] — B [i]  [j];} 

window (170, 200, 650, 480) ;  intensity (B, 250, 250) ;  contour (B, 250, 250, c, nc) 
windowOO;  color(7);  move  (175, 205) ;  printf  ("Max  (%4 . 2f)  ",  0  .  l*Bgrx)  ; 

/ *============  B-z  graphics,  graph  at  r=rl  ==================*/ 

windowO () /window (170, 485, 650,715) ; scale (0, 500, 1000, -500) /color (0) ; 

rect (0,500, 1000,-500) /boldvec (0,0, 1000,  0) /color (4) ; 

boldcatvec (4*i,Bint) ;nrl=(int)  (125.0* (rl+1.5) / 1 . 5 )  ; 

for  (i=0; i<250; ++i)  {Bint-(int)  (  (500 . 0/ (1 . l*Bgrx) ) *B [i]  [nrl] ) ; 

printf ("%d, %d, \n" , 4*i, Bint) ; } 

printf (";  ") /windowO () /color (0) /move (85,  485)  ; 

printf ("%4 . Of", 1 . 1*0 . l*Bgrx) ; 

move (145,  590) /printf ("0") /move (70, 590) /printf ("B (g) «)  ; 
move (75, 700) /printf  ("%4 . Of ", -1 . 1*0 . l*Bgrx)  ; 
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/ *===================  final  ==========================“ - *  / 

windowOO;  color  (0); 

bmn=bmx=0 .0;  for  (i=0;i<250;++i)  {if  (bmn>=B [i] [nrl] )  bmn=B [i] [nrl] ; 

if  (bmx<B [i] [nrl] )  bmx=B[i] [nrl] ; } 

move (675, 310) ;  printf ("an=  %2.2f",an); 

move (675,340) ;  printf ("aT=  %2.2f",aT); 

move  (675,  370)  ;  printf  ("T  bg=  %g  eV",Tbg); 

move (675, 400) ;  printf ("max  T") ; sub ('e' ) ; printf ("=  %g  eV",Tmx); 

move (675, 440) ;  greek (Delta) ;printf ("rn=  %g  cm",l/DZn); 

move (675, 470)  ;  greek (Delta) ;printf ("rT=  %g  cm",l/DZt); 

move (675, 500) ;  printf ("Bmax=  %3.3f  gauss", 0 .l*bmx) ; 

move (675, 530) ;  printf ("Bmin=  %3.3f  gauss", 0 . l*bmn) ; 

move (675, 560) ;  printf ("r=  %g  cm",rl); 

move (675, 590) ;  printf ("dt=  %g  nsec",dt); 

move(675,  620)  ;  printf  ("Z")  ;sub('o' )  ,-printf  ("=  %3.3f  cm",z0); 
mode ("FI"); 

for  (i=0; i<55;++i)  {  color(8+i);  rect  (680+ (int)  (4 . 2*i) , 200,  910, 225) ;  } 
color  (7);  if  (nc!=0)  for  (i=0 ; i<nc; ++i) 

{  it- 680+ (int)  (230.0*  (i+0 . 5) / (f loat ) nc) ;  vector (it, 200, it , 225) ;  } 
mode("F0");  color  (0);  rect  (680, 200,  910, 225) ; 
move (665, 228) ;  printf ("-max") ; 

move (790,228) ;  printf ("B");  move (870, 228) ;  printf ("max") ;  } 


Computer  Simulation  Program  pl2.c. 


B.  PL7.C 


y************************************************* 

====  Self  Generated  Magnetic  Fields  7  ===== 

Simulates  cross  field  density  and  temperature  diffusion 
assuming  fully  ionized  gas  using  Bohm  diffusion. 

PL7.c 

**************************************************/ 

♦include  <stdio.h> 

♦include  <cmath.h> 

♦include  <greeks.c> 

ARRAY 2d  no,  nn,  T; 

float  t , r , nmax, dz , dr , dt , tmax, ngrmx, ngrmn, 

dn,Bf ,A,al,a2,a3,a4,Drn,Tmx,DrT,rl,r2,pi,pf ,Dth,dT 
, Tgrmx, Tgrmn, Ti, Tf , Tbg, nbg, rTe, rne; 

int  i,nz,nr, j,nnt,k,l; 

main ( )  { 

printf ("\n***  Radial  Plasma ' diffusion  ***\n") ; 
printf ( "  tmax,  B,  Drn,  DrT,  Tmx,  Tbg,  nbg,  nmax\n  "); 
scanf  ("%f  %f%f  %f%f  %'f  %f  %f  ",  &tmax,  &Bf ,  &Drn,  &DrT,  &Tmx,  &Tbg,  &nbg,  &nmax) 
dr=5 .0/250.0; dt=0 .001; 

A*=6 . 25e-3/Bf ; Drn*=l/Drn;  DrT=l/DrT; 

for  (nr=0;nr<250;++nr)  { 

r- (2. 5/125.0)* (float) nr; 

T  [5] [nr]=T [5]  [nr]  =  ( (Tmx-Tbg) *exp (-r*r*DrT*DrT) ) +Tbg; 
no [5] [nr] =nn [5] [nr] = ( (nmax-nbg) *exp (-r*r*Drn*Drn) ) +nbg; } 

/ *===========  calc  initial  number  of  particles  ============*/ 

pi=pf=0 .0; 

for  (nr=0;nr<250;++nr)  { 
r=fabs( (2.5/125.0)* (float) nr) ; 
pi=pi+ (6 .28318*r*no[5] [nr] *dr) ; } 

/*================  calc  initial  energy  ====================*/ 

Ti=Tf=0 .0; 

for  (nr=0;nr<250;++nr)  { 
r=fabs( (2.5/125.0)* (float) nr) ; 

Ti=Ti+(6.28318*r*T[5] [nr] *dr) ; } 
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nz=5; 

k= (int) ( (0.832/DrT) *(125.0/2.5) ) ; 
l-(int) ( (0.832/Drn)* (125.0/2.5) ) ; 
for  (t=0 . 0; t<tmax;t=t+dt)  { 

/*==========  Calculate  updated  T  and  n  gradients  ==========*/ 

for  (nr=k;nr<248;++nr)  { 

r=  (2 . 5/125 . 0)  *  (float)  nr;rl=  (2 .5/125 . 0)  *  (float)  (nr+1)  ; 

r2=(2. 5/125.0)* (float) (nr-1); 

al=T [nz] [nr+1] *T [nz] [nr+2] ; 

a2=T [nz] [nr+1] *T [nz] [nr] ; 

a3=T[nz] [nr-1] *T[nz] [nr] ; 

a4=T [nz] [nr-1] *T [nz] [nr-2] ; 

/*==================  Calculate  updated  n-field  =============*/ 

dT= (rl* (al-a2) +r2* (a4-a3) ) / (4 . 0*r*dr*dr)  ; 

T [nz] [nr] =T [nz] [nr] + (A*dT*dt) ; 

} 

/ *============  Calculate  updated  temperature  gradient  ======*/ 

for  (nr=l;nr<248;++nr)  { 

r= (2 . 5/125 . 0) * (float) nr;rl= (2 .5/125 . 0) *  (float) (nr+1) ; 

r2=  (2. 5/125.0)* (float)  (nr-1); 

al=rl*no[nz] [nr+2]; 

a2=rl*no[nz] [nr]; 

a3=r2*no[nz] [nr]; 

a4=r2*no[nz] [nr-2]; 

/*=====  calculate  updated  T-field  and  end  time  loop  ========*/ 

dn= ( (al-a2) +(a4-a3) ) / (4 . 0*r*dr*dr) ; 

Dth=4 .253e-5*pow(T [nz] [nr] ,2.5) /no[nz] [nr] ; 
nn[nz] [nr]=nn[nz] [nr] + (Dth*dn*dt) ; 

) 

for  ( nr=0 ; nr<2 5 0 ; ++nr )  no [nz]  [nr]=nn[nz]  [nr] ; 

} 

/*========  calc  final  n-umber  of  particles  ===================*/ 

for  (nr=0;nr<250;++nr)  { 

r=fabs ( (2.5/125.0)* (float) nr) ; 
pf=pf+ (6 . 28318*r*no [5] [nr] *dr) ; } 

/*==========  calc  final  energy  ==============================*/ 

for  (nr=0 ; nr<250; ++nr)  { 
r=fabs( (2. 5/125.0)* (float) nr) ; 

Tf=Tf+(6.28318*r*T[5] [nr] *dr) ; } 

/*==============  Find  ngr  and  Tgr  max/min  ===================*/ 

ngrmx=ngrmn=Tgrmx=Tgrmn=0 . 0 ; i=5 ; 
for  ( j=0; j<249;++ j)  { 
if  (ngrmx<=no [i] [ j] )  ngrmx=no [i] [ j] ; 
if  (Tgrmx<=T [i] [ j] )  Tgrmx=T [i] [ j] ; 
if  (ngrmn>no [i]  [ j] )  ngrmn=no[i]  [ j]  ; 
if  (Tgrmn>T [i] [ j] )  Tgrmn=T [i] [ j ] ; } 
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init () ;  windowO () ;color_scale ("bluered") /grey_scale ("greyscalel") ; 
bgcol (7) ;  erase  ();  color (0);  mode ("Fl") ;  rect (237, 80, 841, 195) ;  color (2)/ 
rect  (230,75,834,190) ;  modeC'FO");  color  (0);  rect (230, 75, 834, 190) / 

move (270, 80) ;  color  (4);  printf("**  Radial  Cross  Field  Diffusion  **")/ 
color (0) ; 

move (250, 115) ;  printf ("Time=%g  nsec", tmax)  / 

move (250, 140) ;  printf ("max  T" ) / sub (' er ) ; printf ("  =  %g  eV",Tmx); 
move (530, 115) ;  printf ("max  n") ; sub (' o' ) ; printf ("  =  %g  r.u.",nmax); 
move (85,300) ;  printf ("T (eV) ") ; 
move (110, 403) ;  printf  ("0 . 0") / 

move (167,  655) ;  printf ("0") ;  move (380,  655)  ;  printf ("r (cm) ") ; 
move  (605,  655)  ;  printf  ("5.0") ;  windowOQ; 

/*===================  T-r  graphics  ==========================*/ 

windowO ( ) ; window (170,200, 650,420) ; scale (0,1000,1000,0) ; color (0) ; 
rect (0,1000,1000,0) ; color (4) ;i=nnt=0; 
boldcatvec (4*i, nnt) ; 

for  (i=0 ; i<250 ; ++i)  {nnt= (int ) ( (1000 . 0/ (1 . l*Tgrmx) ) *T [5] [i] ) / 
printf ("%d, %d, \n", 4*i,nnt) ; } 
printf ("; ") ; 
mode ("Fl" ) ; 

color (1) /vector (4*k, 1000, 4*k, 0) /mode ("F0") / 
windowO ()  /  color (0) /move (100,200)/ 
printf ("%3. If", l.l*Tgrmx) / 

windowO () / window (170, 425, 6fo,  645) / scale (0, 1000, 1000, 0)  /  color (0) / 
rect (0, 1000,1000, 0) /color (4) /i=nnt=0/ 
boldcatvec (4*i, nnt) / 

for  ( i=0 / i<250/ ++i)  {nnt- (int) ( (1000.0/ (l.l*ngrmx) ) *no[5] [i] ) / 
printf ("%d, %d,\n", 4*i, nnt)/ } 
printf ("/") / 
mode ( "Fl" ) / 

color (1) /vector (4*1, 1000, 4*1,0) /mode ("F0") / 
windowO () /color (0) /move (100, 425)  / 
printf ("%3 . If ",  1 . l*ngrmx) / 

move (108,  628) /printf ( "0 . 0") /move (55, 520) /printf ("n(r.u. ) ") / 


windowO ( ) /  color  ( 0 ) / 

move (675,440) /  printf ("T") /sub('e' ) /printf (" (0)  =  %2.2f  eV",Tgrmx)/ 

move (675, 350) /  printf ("nbg  =  %2.2f  r.u.",nbg)/ 

move (675,380) /  printf ("Tbg  =  %2.2f  eV",Tbg)/ 

move (675,410) /  printf ("n (0)  =  %2.2f  r.u.",ngrmx)/ 

move (675, 515) /  greek (Delta) /printf ("rT") / sub (' o' ) / 

printf ("  =  %3.2f  cm",l/DrT) / 

move (675,480) /  greek(Delta) /printf ("rn") ; sub ('o') / 
printf ("  =  %3.2f  cm",l/Drn)/ 
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move (675, 555) ;  printf("B  =  %g  gauss", Bf); 
move (675, 585) ;  printf("dt  =  %g  nsec",dt); 

move (675,  615) ;  greek (Epsilon) ; sub (' n' )  ;printf ("  =  %3 . 3f ", (pf-pi) /pi) ; 
move (675, 655) ;  greek (Epsilon) ; sub (' T' ) ;printf ("  =  %3 . 3f ", (Tf-Ti) /Ti) ; 
mode ("FI") ; 

for  (i=0; i<250;++i)  if  (no[nz] [i]<=0.368* (nmax+nbg) )  { 
rne= (2 . 5/125.0)*  (float) i; i=250; } 

for  (i=0;i<250;++i)  if  (T  [nz]  [i]<=0 .368* (Tmx+Tbg) )  { 
rTe= (2 . 5/125 . 0) * (float) i; i=250; } 

move (675,320) ;  printf ("T (4cm)  =  %2.2f  eV", T [nz] [247] ) ; 
move (675, 290) ;  printf ("n (4cm)  =  %2.2f  r .u. ",no [nz] [247] ) ; 
move (675, 260) ;  greek (Delta) ;printf ("rT/")  ; 
greek (Delta) ; printf ("rn  =  %1 . 4f ", rTe/rne) ; 

} 


Computer  Simulation  Program  pl7.c. 


Ill 


REFERENCES 


1.  McKee  L.L.,  An  Investigation  of  the  Self-Generated,  Magnetic  Fields 
Associated  with  a  Laser-Produced  Plasma ,  Ph.D.  Dissertation,  Naval 
Postgraduate  School,  1972. 

2.  Bird  R.S.,  The  Pressure  Dependence  of  Spontaneous  Magnetic  Fields  in 
Laser  Produced  Plasmas,  Ph.D.  Dissertation,  Naval  Postgraduate 
School,  1973. 

3.  Allen  C.W.,  Astrophysical  Quantities,  3rd  Edition,  pg.  35-48,  The 
Athlone  Press,  University  of  London,  1976. 

4.  Chen  F.F.,  Introduction  to  Plasma  Physics  and  Controlled  Fusion, 
Vol.I:  Plasma  Physics,  2nd  Edition,  ch.5,  Plenum  Press,  New  York, 
1974. 

5.  Trivelpiece  A.W.,  Principles  of  Plasma  Physics,  pg.83,  McGraw-Hill 
Book  Company,  Inc.,  New  York,  1973. 

6.  Schmidt  G.,  Physics  of  High  Temperature  Plasmas,  2nd  Edition, 
Academic  Press,  New  York,  1979. 

7.  Stamper  J.A.  and  others,  Spontaneous  Magnetic  Fields  in  Laser- 
Produced  Plasmas,  Phys.  Rev.  Letters,  v.26,  no.17,  pg.  1012-1015, 
1971. 

8.  Dean  S.O.,  McLean  E.A.,  Stamper  J.A.,  and  Grien  H.R.,  Phys.  Rev. 
Letters,  v.27,  no.487,  1971. 

9.  Brooks  K.M.,  An  Investigation  of  Early  Disturbances  Found  in 
Association  with  Laser-Produced  Plasmas,  M.S.  Thesis,  Naval 
Postgraduate  School,  1973. 


112 


10.  Schwirzke  F.,  Bird  R.S.,  McKee  L.L.,  Cooper  A.W.,  Pressure 
Dependence  of  Self-Generated  Magnetic  Fields  in  Laser-Produced 
Plasmas,  Phys.  Rev.  A,  v.7,  no.  4,  1973. 

11.  Schwirzke,  F.,  Measurements  of  Spontaneous  Magnetic  Fields  in 
Laser-Produced  Plasmas,  Laser  Interactions  and  Related  Plasma 
Phenomena,  Vol.3,  Plenum  Press,  1974. 

12.  Birdsall  C.K.,  Longdon  A.B.,  Plasma  Physics  via  Computer  Simulation, 
ch.4,  Adam  Hilger  Press,  Bristol,  1991. 

13.  Intrilligator  D.S.,  Miller  W.D.,  Journal  of  Geophysical  Research,  v.87, 
pg.  4354,  1982. 

14.  Glasstone  Samuel,  Dolan  P.J.,  The  Effects  of  Nuclear  Weapons,  3rd 
Edition,  ch.3&ll,  US  Government  Printing  Office,  Washington  D.C., 
1977. 

15.  Book  David  L.,  NRL  Plasma  Formulary,  pg.28-32,  Naval  Research 
Laboratory,  Washington  D.C.,  1980. 


113 


INITIAL  DISTRIBUTION  LIST 

1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5101 

3.  Professor  Fred  Schwirzke,  Code  PH/SW  3 

Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4.  Professor  William  Colson,  Code  PH/CW  1 

Chairman,  Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Professor  A.W.  Cooper,  Code  PH/CR  1 

Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  LT  John  P.  Carter  3 

613  Colonial  Manor  Road 

N.  Huntingdon,  Pennsylvania  15642 


114 


